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Section  1 


INTRODUCTION 


We  have  undertaken  the  development  of  parameterization  programs 
involving  the  transfer  of  thermal  IR  and  solar  fluxes  in  clear  and 
various  cloudy  atmospheres.  The  primary  objective  of  this  research  work 
is  to  construct  a  unified  and  coherent  radiation  package  that  is  computa¬ 
tionally  economic  for  incorporation  into  the  global  models  developed  at 
the  Air  Force  Geophysics  Laboratory  (AFGL).  In  particular,  our  major 
emphasis  is  on  the  development  of  accurate  and  stable  radiation  programs 
for  atmospheres  containing  high,  middle,  and  low  clouds  and  combinations 
of  these  clouds.  The  other  objective  is  to  develop  cloud  formation 
schemes  for  the  GCM  in  connection  with  radiation  programs.  Our  effort  in 
the  design  of  radiation  and  cloud  programs  is  closely  in  line  with  the 
overall  AFGL  effort  in  developing  a  global  model  for  short  and  medium 
range  weather  prediction  and  cloud  forecast  research. 

In  Section  2  of  this  report,  we  describe  the  physical  concepts  and 
basic  parameterization  equations  developed  for  the  computation  of  radia¬ 
tive  fluxes  in  clear  atmospheres.  In  Section  3  we  first  present  the 
approaches  that  are  used  in  predicting  and  prescribing  cloud  properties. 
In  particular,  the  strategy  for  strapping  cloud  layers  into  cloud  decks 
and  the  method  for  determining  total  cloud  cover  are  discussed  in  some 
detail.  We  then  present  the  manner  in  which  the  transfer  of  IR  and  solar 
fluxes  in  multilayered  cloudy  atmospheres  is  parameterized.  In  Section 
comparisons  of  IR  cooling  and  solar  heating  rate  profiles  for  various 
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cloud  configurations  computed  from  the  present  parameterizations  with 
those  from  a  more  exact  doubling/adding  method,  also  developed  at  the 
University  of  Utah,  are  carried  out.  The  present  radiation  and  cloud 
programs  are  successfully  incorporated  in  the  AFGL  12-  and  1 8-layer 
models.  Simulations  are  carried  out  for  cloud  and  radiation  fields  using 
the  AFGL  12-layer  model  and  initial  humidity  and  temperature  fields 
corresponding  to  12Z,  17  June  1979.  Results  of  zonally  averaged  cloud 
cover,  radiation  budgets  at  the  top  of  the  atmosphere  and  surface,  as 
well  as  atmospheric  cooling  profiles  generated  from  the  model  are 
presented  and  discussed  in  Section  5.  Finally,  a  summary  is  given  in 
Section  6. 


Section  2 


RADIATIVE  TRANSFER  IN  CLEAR  ATMOSPHERES 


2. 1  Infrared  Radiative  Transfer 

The  upward  and  downward  fluxes  at  a  given  level  In  a  clear  atmo¬ 
sphere  In  the  path  length  coordinate  may  be  expressed  In  the  forms  (Liou, 
1980) 


F* (u)  =  oT^  [l  -  ef(u,T  )]  +  f  ot\u')  K(u-u')  du' 

S  3  }Q 

*  fU  U 

F  (u)  =  oT  (u' )  K(u'-u)  du'  , 

•  i-\ 


(2.1  ) 

(2.2) 


where  Tg  is  the  surface  temperature,  u  the  path  length  for  the  absorbing 
gas,  and  u1  the  total  path  length.  The  flux  expressions  In  the  u- 
coordlnate  can  be  easily  transferred  to  the  height  or  pressure  coordi¬ 
nate.  The  broadband  flux  emlsslvlty  Is  defined  by 


CD 

ef(u,T)  -  irB^T)  h  '  /(«)]  dv/(o/)  , 


(2.3) 


4> 

where  Tk  Is  the  flux  transmittance  at  the  spectral  wave  number  v  and  ttBn1 
v  v 

the  Planck  flux.  The  kernel  function  K(u-u' )  Is  related  to  the  broadband 
flux  emlsslvlty  as  follows: 


K(u-U')  =  -  i  Ef  [ |u-u’ | ,T(u’ ) ] 


(2.U) 


To  apply  the  broadband  flux  emlsslvlty  for  flux  and  heating  rate 
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calculations,  the  IR  spectrum  is  divided  into  five  bands,  including  (1) 
H20  rotational  band,  (2)  H20  continuum  band,  (3)  H20  vibrational-rota¬ 
tional  band,  (4)  0^  rotational-vibrational  band,  and  (5)  C02  rotational 
vibrational  band.  The  overlap  between  bands  (1)  and  (5)  is  denoted  as 
band  6.  Thus,  the  broadband  flux  emisslvity  may  be  written  m  the 
form 


e1 (u,T) 


e[(VT)  -  ^VVT) 


_  — 


(2.5) 


with  the  individual  broadband  flux  emissivity  defined  by 


,oo 

ef  =  irB-(T)  [l  -  /  .  (u)]  dv/U^J  , 
i  J0  v  v,i 

In  Eq.  (2.5),  u1  =  u2  «  u^  *  uM,  the  H20  path  length,  uij  =  uQ,  the  0^ 
path  length,  u5  =  uc,  the  C02  path  length,  and  ut  denotes  the  pressure 
and  temperature  corrected  path  length  for  the  respective  absorbing 

gaseous  component.  The  correction  due  to  overlap  is  approximated  by 

f  f/~  \  , 

e6  “  -  e5(uc^  • 

The  broadband  emissivlties  for  the  four  bands  may  be  expressed  in 
terms  of  polynomial  functions  in  the  form  (Liou  and  Ou,  1981) 


e 


[(ut,T) 


3 

exp  (  l 
n=0 


'ni 


(2.6) 


where 


U1  =  log10  “l  "  ai^/bi 


a.  =  log,,  (u.  •  u.  ,  j 

i  10  i,max  i,min' 
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The  coefficients  cnj_  were  obtained  by  the  least-square  method.  The 
quantity  is  expressed  as  a  linear  function  of  logi0  u^ ,  so  that  -1  < 

Uj  <  1.  Errors  are  smaller  for  the  intermediate  u^  values  than  for  - 
-1  or  u,  -  1 ,  due  to  the  characteristic  of  the  least-square  method 
(Hemming,  1973).  For  water  vapor  rotational  and  virbrational-rotatlonal 
bands,  the  broadband  flux  emlsslvltles  are  derived  from  the  random  band 
parameters  assuming  the  Lorentz  line  shapes  obtained  by  Rodgers  and 
Walshaw  (1966)  in  which  an  empirical  temperature  correction  for  the 
rotational  band  was  included.  For  the  water  vapor  continuum  absorption 
in  the  window  region,  the  empirical  formula  developed  by  Roberts  et  al. 
(1976)  was  used.  For  ozone,  random  model  parameters  for  the  9.6  pm  band 
derived  by  Goldman  and  Kyle  (1968)  were  utilized  to  derive  the  0^ 
broadband  emissltivlty .  The  prescription  of  Uj  max  and  u^  mln  was 
according  to  the  range  for  realistic  atmospheres.  For  water  vapor, 

Ui.max  =  10  g  cm-2  and  u1>mln  -  10-^  g  cm"2.  Table  1  lists  the  coeffi¬ 
cients  of  the  broadband  emlsslvity  equation  cn±  for  the  three  water  vapor 
(1  *  1,2,3)  and  ozone  bands  (i  -  5). 

The  broadband  emlsslvity  for  carbon  dioxide  is  computed  using  the 
parameterization  scheme  developed  by  Ou  and  Liou  (1983).  In  that  paper 
the  C02  broadband  emlsslvity,  which  is  a  function  of  the  pressure  and 
temperature  corrected  path  length  uQ  and  the  atmospheric  temperature,  is 
derived  from  line-by-line  transmittance  data  presented  by  Fels  and 
Schwarzkopf  (1981).  From  Fig.  1  in  the  paper  by  Ou  and  Liou,  the  broad¬ 
band  flux  emlsslvity  for  the  C02  15  pm  band  can  be  fitted  into  two 
functional  forms: 


A."'.'.''  . v; 


Table  1.  Coefficients  of  broadband  emissivity  values  for  water  vapor 


and  ozone  bands. 


1 

c0i 

c1i 

c2i 

c3i 

T=203°K  1 

2 

3 

-.16378+001 

-.98069+001 

-.583119+001 

-.67785+001 

.31523+001 

.91563+001 

.54245+001 

.42930+001 

-.28231+001 

-.54110+000 

-.34050+001 

-.95681+000 

.10998+001 

-.27859+000 

.85649-002 

-.55536+000 

T=233°K  1 

2 

3 

l* 

-.18865+001 

-.95159+001 

-.50133+001 

-.63805+001 

.32475+001 

.91534+001 

.53027+001 

.42918+001 

-.27705+001 

-.52756+000 

-.34506+001 

-.95672+000 

.10903+001 

-.26671+000 

.10991-001 

-.55860+000 

T=263°K 

1 

2 

3 

4 

-.21226+001 

-.93430+001 

-.44264+001 

-.61317+001 

.33251+001 

.91512+001 

.53278+001 

.42910+001 

-.27277+001 

-.51723+000 

-.34806+001 

-.95628+000 

.10832+001 

-.25749+000 

.79503-002 

-.56004+000 

T=293°K 

1 

-.23454+001 

.33888+001 

-.26927+001 

.10777+001 

2 

-.92462+001 

.91495+001 

-.50918+000 

-.25017+000 

3 

-.39991+001 

.53068+001 

-.35005+001 

.42209-002 

4 

-.59797+001 

.42904+001 

-.95611+000 

-.56109+000 

exp 

3 

(  I  a 

0  n 

u^n)  ,  for  u  ^  10  14 

g 

-2 

cm 

_ 

n=0 

^(u5.T)  =  < 

tb0  *  b1 

u£  J  ,  for  u  <  10 

-2  >  <2-7 

exp 

g 

cm 

where 

u ' 

’  -  (2 

1,eto  “ 

+  7 . 69897  ) /6 . 301 03  . 

The  coefficient  values 

are  aQ 

=  -4.00893  +  f (T) ,  a1 

s 

4.39828,  a p  = 

-3.07709,  a3  -  0.94529,  b0  -  -**.0036  +  f(T),  and  bj  =  5.13453,  where 


f(T)  =  in  {h(T)  g(T)/[h(T0)  g(T0)]}  . 


The  quantities  h(T)  and  g(T)  are  given  by 

1  -  T(u  ,t) 

h(T)  =  -  =  1  +  AAT  (1  -  BAT )  ,  (2.8) 

1  -  T(u5.T0) 

fV2  2j 

g(T)  *  irB  (T)  dv  /  o T  ,  (2.9) 

V1 

where  T(u5,t)  is  the  transmittance,  AT  =  T  -  T0,  A  =  1  .833x10“\  B  *= 
1.364x10“2,  Tq  =  250°K,  v-j  =  500  cm-1,  and  v2  -  850  cm-1.  Only  the 
coefficients  aQ  and  bQ  are  temperature-dependent  since,  by  using  the 
definition  of  Hu^.T)  (Fels  and  Schwarzkopf,  1981),  z~  can  be  expressed 
by 


ef(u5,T)  =  h(T)  g(T)  [l  -  2*(u5,T0)] 


(2.10) 


In  this  manner,  the  dependence  of  the  temperature  T  and  path  length  uc  on 
the  emisslvity  has  been  separated. 

The  pressure  and  temperature  adjustments  on  the  path  length  for  each 
band  are  different.  For  the  three  water  vapor  bands,  pressure- 
temperature  scaling  methods  proposed  by  Llou  and  Ou  (1981)  are  used. 

These  are 


U1 


P(Z,)  exp  |A'  [T(z')  -  T  ]  +  B'  [ T(z ' )  -  T  ]2}  p(z')  dz*  , 
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exp 
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(2.11  ) 


1800  r  ,  I 

T-n^ry  [T(z’)  -  Tb]  S  pw(z-)  dz'  ,  (2.12) 
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\h\ 


I 

•«! 
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-  r p(z,)  r?o  T/2 

u3  ’  J0  PQ  [W7J  P* 


( z ' )  dz'  . 


(2.13) 


where  A»  =  10.6278x10"3  K"1 ,  B'  -  -44.6152x10“°  K-2,  Ta  -  260°K,  Tb  - 
296°K,  Tc  ■=  300°K,  P0  -  1013  mb,  Pw  and  P0w  are  vapor  pressures  at 
temperatures  T  and  Tb,  respectively,  and  pw  is  the  water  vapor  density. 

The  pressure  and  temperature  corrected  path  lengths  for  COp  have  the 
form  (Ou  and  Llou,  1983) 


uc  =  2  c  /  {[l  +  M(c2/u2  +c/u)]1/2  -  l}  , 
o 


_h  _0 

where  the  constant  c  -  3*7551x10  g  cm  and 


(2.U) 


[Z  P(z' )  T0 

U5  '  J0  F0  FrT| 


p  (z* )  dz' 


with  pc  the  carbon  dioxide  density  and  Tq  »  273°K.  In  Eq.  (2.14)  a 
combination  of  u  and  u  is  used,  since  we  were  unable  to  find  a  single 
pressure-scaled  path  length  that  could  provide  accurate  flux  emissivity 
values.  Finally,  no  adjustment  has  been  made  for  the  ozone  path  length 
because  of  the  uncertainty  In  applying  pressure  and  temperature 
corrections. 


2.2  Solar  Radiative  Transfer 

In  a  clear  atmosphere,  Rayleigh  scattering  due  to  air  molecules  is 
insignificant  above  a  height  zy  ,  about  10  km.  The  upward  solar  flux 
above  this  height,  on  the  other  hand,  is  contributed  from  upward  Rayleigh 
reflection  plus  a  diffuse  reflection  component  that  results  from  multiple 
reflections  between  the  Rayleigh  layer  and  underlying  surface.  Thus, 


TO 


the  downward  and  upward  solar  fluxes  at  level  z  above  z1  may  be  expressed 
by 

F+(z)  -  p0  SQ  [l  -  A+(zt-z)]  ,  (2.15) 

F+(z)  =  F*(Zl)  (rg  +  G  ta)  [l  -  r(z-Zl)]  ,  (2.16) 

where  zT  is  the  height  at  the  top  of  the  atmosphere,  the  solar  constant 
S0  Is  taken  to  be  1360  W  m~2,  and  the  cosine  of  the  solar  zenith  angle  is 
given  by  (Liou,  1980) 

Uq  =  sine  slnfi  ♦  cose  cos6  cosA  ,  (2.17) 

where  6  is  the  latitude,  6  the  inclination  angle,  which  is  a  function  of 
the  day  of  the  year,  and  A  the  hour  angle.  The  flux  reflection  ra  due  to 
Rayleigh  scattering  is  taken  to  be  ra  =  0. 28/(1  +  6. *13  p0)  (Lacis  and 
Hansen,  197*0.  The  flux  transmission  is  given  by  ta  =  1  -  ra.  ra  and  ta 
are  the  corresponding  diffuse  (or  global)  values  that  are  obtained  from 
integration  over  the  positive  y0  values.  The  direct  downward  gaseous 
broadband  solar  absorptivity  A^(z)  will  be  defined  later  and  its  diffuse 
value  A+(z)  »  A+  [n  u(z)]  where  the  diffusivity  factor  n  =  1.66  for  H2O 
and  1.9  for  0^.  The  higher  diffusivity  factor  for  ozone  is  chosen  to 

account  for  Rayleigh  scattering  effects  on  the  ozone  absorption.  Based 

on  a  simple  adding  procedure,  the  multiple  reflections  involving  the 
Rayleigh  scattering  layer  and  the  surface  having  a  Lambertian  albedo  of 
rg  are  described  by  a  nondimensional  upward  diffuse  component  given 


In  Eq.  (2.16),  the  ra  term  represents  the  reflection  of  the  direct  down¬ 
ward  solar  flux,  while  G  ta  denotes  the  diffuse  reflection  component. 
Likewise,  downward  and  upward  fluxes  at  heights  below  z1  may  be  expressed 
by 

F*(z)  -  uQ  SQ  [1  -  A+(zt-z)]  ,  (2.19) 

F*(z)  -  F+(0)  r  [l  -  A * ( z )  ]  .  (2.20) 

s 

It  should  be  noted  that  Rayleigh  scattering  below  z-,  will  reduce  both 
upward  and  downward  fluxes  by  about  the  same  amount  since  its  overlap 
with  gaseous  absorption  in  the  near  infrared  may  be  neglected.  Thus,  the 
net  flux  and  hence  the  heating  rate  may  be  evaluated  without  the  incorpo¬ 
ration  of  Rayleigh  scattering. 

Three  absorbers  are  considered  in  the  transfer  of  solar  radiation, 
i.e.,  water  vapor,  carbon  dioxide,  and  ozone  with  the  former  two  gases 
responsible  for  absorption  in  the  near  infrared  region  and  the  latter 
primarily  in  the  UV  region.  The  overlap  of  H20  and  0^  absorption  is 
insignificant  so  that  the  broadband  solar  absorptivity  may  be  expressed 
by 

A(z)  =  l  aJu^Uq)  Ti  ♦  I  A1(uw/u0)  ♦  e  A5(u/m0)  ,  (2.21) 

where  uQ,  uw,  and  uc  represent,  respectively,  path  lengths  for  0^,  H20, 
and  C02,  and  f^  are  the  fractional  solar  fluxes  associated  with  the 
individual  band  absorptl vities  Aj.  The  H20-C02  overlap  correction  is 
given  by  c  »  a  +  b  A^(uw)  where  a  and  b  are  empirical  constants 
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derived  from  the  measured  data  presented  by  Howard  et  al.  (1956);  a 
0.75  and  b  =  -0.0015  for  A5(uw)  A\>5  l  300  cm”1,  and  a  -  1  and  b  - 
-0.00233  for  A5(uw)  <  300  cm”1. 


The  Individual  absorptivity  for  ozone  is  given  by 

AJuo/po)  “  1  ‘  exp  t-ki  Uo/uo)  ’ 


(2.22) 


where  k^  denote  the  ozone  absorption  coefficients  taken  from  the  data 
obtained  by  Inn  and  Tanaka  (1953)  and  Vigroux  (1953)*  and  uQ  is  in  units 
of  cm-atm.  Table  2  lists  the  values  of  the  ozone  absorption  coefficients 
and  fractional  solar  fluxes. 

For  water  vapor  and  carbon  dioxide,  the  empirical  band  equation 
derived  by  Llou  and  Sasamori  (1975)  is  used.  It  has  the  following 


analytic  form 


si<u/“o!  -  *7  ci  *  Di 


K./D.  -C./D. 

l°g,0  —  P  *  10 


,  (2.23) 


where  Av^  are  the  band  widths,  C^,  D^,  and  are  empirical  constants. 
The  pressure  P  in  units  of  mm  Hg  is  weighted  over  the  path  length  given 


-  f  P(u 
J0 


)  du/u  , 


_  p 

where  u  can  be  either  uw  or  uc  and  is  in  units  of  g  cm  .  The  values  of 
Ct,  Dlt  and  K1  for  six  H20  and  the  2.7  pm  C02  bands  are  listed  in  Table 
3,  along  with  the  fractional  solar  fluxes. 
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Table  2.  Ozone  absorption  coefficients  and  fractional  solar  fluxes. 


(-14H°C) 


3  0.20-0.21 

1 

9.8 

1  -  244x1  0~JJ 

|  0.21-0.22 

2 

27 

2.97x10“^  ' 

1  0.22-0.23 

3 

75 

44.59x10*2  ! 

<  0.23-0.2^4 

44 

1644 

4J.59X10"4 

J  0.244-0.25 

5 

2544 

5.144x10*2 

2  0.25-0.26 

6 

290 

7.55x10"^ 

ft  0.26-0.27 

7 

24»1 

1.35x10*3 

0.27-0.28 

8 

14)5 

1 .59x10"3 

0.28-0.30 

9 

33.7 

6.446x10*3  ! 

0.30-0.32 

10 

2.8 

I.OIxlO*2 

0.32-0. 3^ 

11 

0.16 

1 .50x10*2 

o.3n-o.35 

12 

0.0144 

7.95x10  3 

0.445-0.50 

13 

0.0107 

7.446x10‘2 

r  0.50-0.55 

in 

0.055 

6.78x1 0~2 

0.55-0.60 

15 

0.11 

6. 30x1 0*2  ^ 

0.60-0.65 

16 

0.09 

5.87x1 O*2  j 

$  0.65-0.70 

17 

0.038 

5.33x10"2  $ 

I  0.70-0.80 

18 

0.015 

9 . 1 4Jx  1 0*2  | 

Table  3.  Empirical  constants  for  H20  and  C02  bands  and  fractional  solar 


fluxes. 


Section  3 


RADIATIVE  TRANSFER  IN  MULTILAYERED 
CLOUDY  ATMOSPHERES 

In  this  section,  we  shall  first  present  the  methods  for  predicting 
or  prescribing  the  required  cloud  properties  as  input  for  radiation 
calculations.  We  then  describe  the  parameterizations  developed  for  com¬ 
puting  both  IR  and  solar  flux  and  heating  rate  profiles  in  atmospheres 
containing  multiple  stratiform  cloud  layers. 

3. 1  Cloud  Property  Prediction/Prescription 

Radiation  computations  require  the  cloud  top  and  base  heights,  cloud 
cover,  and  cloud  water  content  as  input  parameters.  In  the  present 
radiation  model,  the  vertical  structure  is  the  same  as  that  used  in  the 
general  circulation  model  (GCM),  l.e.,  both  the  radiation  model  and  GCM 
have  been  modified  to  accomodate  variable  resolution.  Due  to  the  method 
used  in  predicting  the  layer  cloud  cover,  each  model  layer  has  a  possi¬ 
bility  for  cloud  formation.  We  then  "strap"  two  or  more  contiguous 
layers  into  cloud  "decks."  At  present,  the  model  allows  at  most  three 
decks.  In  each  deck,  the  cloud  is  allowed  to  occupy  a  fraction  of  the 
whole  horizontal  grid  space.  In  the  vertical,  however,  the  cloud  is 
assumed  to  fill  the  whole  deck  domain.  Thus,  the  cloud  top  is  located  at 
the  top  of  the  highest  model  layer  in  the  deck,  and  the  cloud  base  is  at 
the  bottom  of  the  lowest  model  layer  in  the  deck.  This  is  illustrated  in 
Fig.  1.  In  the  radiation  computation,  we  assume  that  all  cloud  decks 


Fig . 1  (a)  Partial  cloudiness  is  formed  in  each  model  layer, 

clouds  are  strapped  into  three  decks  of  clouds. 


occupy  the  same  fraction  In  the  horizontal  grid  space,  and  we  call  this 
fraction  the  "total  cloud  cover."  This  assumption  enables  us  to  save 
computation  time,  as  the  number  of  vertical  integrations  for  radiation 
computations  is  reduced  to  two  in  a  partially  cloudy  column.  In  general, 
however,  the  total  cloud  cover  may  be  determined  by  certain  overlap 
assumptions.  In  the  remaining  part  of  this  subsection,  we  shall  describe 
the  cloud  cover  prediction  method,  the  strapping  of  cloud  layers,  and  the 
assumptions  used  to  determine  the  total  cloud  cover.  The  cloud  water 
content,  on  the  other  hand,  is  prescribed  according  to  typical  observed 
values  (e.g.,  Mason,  1971). 

3.1.1  Layer  Cloud  Cover 

Most  available  cloud-cover  prediction  schemes  are  diagnostic.  They 
all  have  similar  parameterized  forms  (Geleyn,  1981;  Slingo  and  Ritter, 
198H,  etc.): 


But  they  differ  in  the  prescription  of  the  threshold  relative  humidity 

V 

In  our  present  program,  we  used  Geleyn’s  specification  of  hc,  i.e., 
h  ■  1  -  ao  (1-o)  [1  +  6 ( 0- 1/2)3  ,  (3*2) 

where  a  =  2,  8  *  /3,  0  =  p/p#,  and  p#  is  the  surface  pressure.  In  the 
AFGL  12-layer  GCM,  clouds  are  formed  in  layers  3~9.  The  top  two  layers 
and  the  bottom  three  layers  are  specified  as  being  cloudless. 


3.1.2  Cloud  Strapping  Scheme 

Based  on  surface  observations,  there  are  generally  less  than  three 
cloud  decks  In  the  vertical.  Accordingly,  It  appears  appropriate  to 
strap  the  model-generated  multilayer  cloudiness  into  at  most  three  cloud 
decks.  Since  in  most  GCMs,  the  vertical  structure  for  both  radiation  and 
dynamic  programs  is  fixed,  it  would  be  extremely  difficult  to  allow 
partial  cloudiness  in  the  vertical  direction.  Hence,  as  we  strap  the 
clouds  into  the  cloud  deck,  it  is  assumed  to  fill  the  respective  deck 
domain  vertically.  For  the  AFGL  12-layer  GCM,  the  high,  middle,  and  low 
cloud  decks  fill,  respectively,  layers  (3,  *0 ,  layers  (5,  6),  and  layers 
(7,  8,  9).  In  this  manner,  the  distinction  of  cloud  decks  is  somewhat 
arbitrary.  Nevertheless,  the  boundary  between  high  and  middle  clouds 
roughly  corresponds  to  the  -20°C  isotherm,  which  has  been  used  as  a 
boundary  between  these  clouds  by  other  researchers  (e.g.,  Koenig  et  al . , 
1987). 

3.1.3  Total  Cloud  Cover 

The  t,otal  cloud  cover  for  each  horizontal  grid  space  is  determined 
in  two  steps.  The  first  step  involves  calculations  of  the  cloud  cover 
for  each  deck.  Then,  the  total  cloud  cover  may  be  evaluated  by  using 
certain  overlap  assumptions.  There  are  two  schemes  that  car  be  used  to 
deterrr.ire  the  cloud  deck  cover, 
i)  Average  Scheme: 

The  cloud  cover  for  each  deck  of  clouds  is  the  average  of  the 
respective  cloud  cover  of  model  layers.  For  example,  using  the  AFGL  12- 
layer  model  structure,  we  have 


(3.3) 


\  “  7  K  +  ' 


nm  "  7  (n5  *  ’ 


nil  "  T  K  +  n8  +  n9^  ’ 


where  nh,  nm,  and  n^  denote  cloud  covers  for  high,  middle,  and  low 


clouds. 


ii)  Maximum  Scheme: 


The  cloud  cover  for  each  deck  of  cloud  is  the  maximum  of  the 
respective  cloud  cover  of  model  layers.  Thus,  for  the  AFGL  12-layer 
model,  we  have 


nh  -  max  (n^.hjj)  , 
nm  -  max  (ryn6)  , 

*  max  (ny,ng,hg)  . 


(3.H) 


In  the  present  model,  we  have  used  the  maximum  scheme. 

To  determine  the  total  cloud  cover,  we  must  consider  the  effect  of 
overlap  between  cloud  decks.  Several  assumptions  can  be  made  on  the 
overlap  configurations. 

1)  Maximum  Overlap 

In  this  case,  the  total  cloud  cover  is  given  by 


n  -  max  (ri!(*r\n,nh) 


(3.5) 


11 )  Random  Overlap 

This  configuration  assumes  that  clouds  overlap  each  other  In  a 
statistically  random  manner.  Thus,  the  total  cloud  cover  may  be 
expressed  by 
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n  -  1  -  (i  -  n^]  (i  -  nm)  (i  -  nh)  . 


(3.6) 


Hi)  Minimum  Overlap 


1 

t 


m 
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This  configuration  assumes  that  clouds  overlap  each  other  in  a 
minimum  fashion,  so  that  there  are  only  one-layer  cloud  configurations, 

i  .e. , 


n  =  min  ^  *  hm  ♦  nh)  .  (3.7) 

It  appears  reasonable  to  assume  that  the  cloud  configuration  ir 
convectively  active  and  overcast  regions  would  tend  to  be  of  maximum 
overlap.  In  the  region  of  small  cloudiness,  on  the  other  hand,  the 
cloud  configuration  would  tend  to  be  of  minimum  overlap. 

In  the  present  program,  we  use  the  maximum  overlap  assumption.  In 
each  partially  cloudy  grid  space,  the  total  cloud  fraction  n  =  max 
(n^.nm.ih)  and  the  clear  fraction  is  1  -  n.  Although  each  cloud  deck  has 
a  different  cloud  cover,  for  the  sake  of  saving  computer  time,  we  set 
nfc  =  nm  “  nh  =  1,  whenever  na,  nn,  nh  >  0. 

3.2  Parameterization  of  Radiative  Transfer  in  Cloudy  Atmospheres 
3.2.1  Cloud  Configurations 

Based  on  the  scheme  described  previously,  there  are  seven  cloud 
configurations.  These  include  three  single-decks,  three  two-decks,  ard 
one  three-deck.  Using  maximum  overlap,  we  first  determine  the  configu¬ 
ration  that  fits  the  grid  point  of  interest.  Then  the  cloud  coven  of 
each  deck  is  assumed  to  be  the  same  as  the  total  cloud  cover.  If  the 
total  cloud  cover  is  less  than  0.1,  the  grid  is  considered  to  be  clear. 
If,  however,  the  total  cloud  cover  is  greater  than  0.9,  the  grid  is 
courted  as  completely  overcast.  For  the  clear  and  overcast  cases,  only 


$ 
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one  vertical  column  radiation  calculation  is  necessary.  For  partially 
cloudy  conditions,  two  calculations  are  required. 


3.2.2  IR  Flux  Exchange 


We  consider  a  three-deck  cloud  system,  as  shown  in  Fig.  2.  The  high 


cloud  deck  is  considered  to  be  nonblack.  Let  the  broadband  emisslvlty, 


transmissivity,  and  reflectivity  for  the  high  cloud  deck  be  denoted  by 


ec ,  tc,  and  rc ,  respectively.  In  a  manner  similar  to  that  documented  in 


Llou  et  al .  (1984),  the  flux  components  for  each  clear  region  separated 


by  three  decks  of  clouds  may  be  derived.  Above  the  high  cloud  we  have 


*  fZ  4 

F  (z)  »  oT  (z')  K(z’-z)  dz'  , 

>7 


(3.8) 


F+(z)  -  {[r°  F*(zt1)  ♦  tC  F*(zbl)  -  (l-eCJ  oT^  ] }  (l  -  /  ( z-zfc1  ,TU  )  ] 


oTt1  [l  -  e  (z-z  ,T  )]  +  oT  (z’)  K(z-z')  dz’  ,  (3.9) 

*  •y 


where  Zj  denotes  the  top  height  in  the  model,  Tt1  the  high  cloud  top 
temperature,  and  zb^  and  z^  the  high  cloud  top  and  base  heights,  respec¬ 


tively.  The  expression  for  the  first  term  on  the  right-hand  side  of  Eq. 


(3.9)  is  explained  in  Appendix  A.  Below  the  high  cloud  and  above  the 


middle  cloud  we  have 


F*(z)  -  lrC  F*Um)  .  tC  F4Ut,)  -  (l-E°)  dj,]  (1  -  /(ZtJ,-z,Tb1)l 


*  oTbl  h  ■  c‘Ub,-z,Tb1)J  *  [  oT#(z*)  K(z'-z)  dz  ,  (3.10) 
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IR  Transfer 


F  (z)  =  oTt2  [l  -  e  U-zt2,Tt2] J  +  j  oT  (z')  K(z-z’)  dz'  , 


(3.11  ) 


where  Tb1  is  the  high  cloud  base  temperature,  Tt2  the  middle  cloud  top 


temperature,  and  zb2  the  middle  cloud  top  height.  Again,  the  first  term 


on  the  right-hand  side  of  Eq.  (3.10)  is  explained  in  Appendix  A.  Both 


middle  and  low  clouds  are  considered  to  be  blackbodies.  Below  the 


middle  cloud  and  above  the  low  cloud  we  have 


+  ij  r 

F  (z)  =  oTb2  [1 


fZ  U 

(zb2~z,Tb2)]  +  oT  (z')  K(z-z')  dz'  , 

~  * 


(3.12) 


F+  (z )  =  oT^  [l  -  ef(z-zt3,Tt3)]  ♦  aAz')  K(z-z' )  dz’  , 


(3.13) 


where  Tb2  is  the  middle  cloud  base  temperature,  the  low  cloud  top 
temperature,  zb2  the  middle  cloud  base  height,  and  zfc3  the  low  cloud  top 


height.  Below  the  low  cloud  and  above  the  surface  we  have 


F+(z)  =  oTb3  [l  -  ef'(zb3-z,Tb3)  j  +  oT^z')  K(z-z')  dz'  , 


(3.14) 


z 

Ff(z)  =  oT^  [l  -  e1 (z,T  )]  +  oT^tz')  K(z-z' )  dz'  ,  (3.15) 

s  3  n 


where  Tb^  is  the  low  cloud  base  temperature,  Ts  the  surface  temperature, 


and  zb3  the  low  cloud  base  height, 


When  there  is  only  one  cloud  deck  generated  in  the  model,  it  can  be 


either  high  (nonblack),  middle  (black),  or  low  (black).  If  the  cloud  is 


nonblack,  we  set  zt2  -  o,  and  Tt2  =  Ts,  and  use  Eqs.  (3.8M3.11)  only. 

If  the  cloud  is  black,  we  use  Eqs.  (3.8)-(3.11)  with  rc  -  0,  ec  =  1  , 
zt2  =  0,  and  T^2  =  Ts.  When  there  ane  two  cloud  decks  generated  In  the 
model,  the  upper  deck  can  be  either  high  (nonblack)  or  black.  In  the 
case  of  a  high  cloud  on  top  of  a  black  cloud,  we  set  zt3  =  0  and  Tt3  = 

Ts,  and  use  Eqs.  (3-8)— (3-11 )-  In  the  case  of  two  decks  of  black  clouds, 
we  use  Eqs.  ( 3. 8)- (3- 1 3)  with  rc  =  0,  ec  «=  1  ,  zt3  =  0,  and  Tt3  =  Tg. 
Furthermore,  if  all  three  decks  of  clouds  are  present  in  the  model,  the 
entire  set  of  Eqs.  (3-8)— (3-15)  is  used  in  the  IR  flux  calculation. 

The  flow  diagram  illustrating  the  strategy  for  distinguishing 
between  clear,  one-,  two-,  and  three-deck  cloudy  atmospheres  is  shown  in 
Fig.  3.  The  integral  quantitites  LC(1),  LC(3),  and  LC(5)  denote  the 
cloud  top  level  indices  corresponding  to  cloud  decks  one,  two,  and  three, 
respectively.  On  the  other  hand,  LC(2)  denotes  the  cloud  base  level 
index  for  cloud  deck  one,  and  LCTOP  and  LCBOT  represent  the  high  cloud 
top  and  base  level  index,  respectively.  We  set  LC(1),  LC(3),  and  LC(5) 
to  zero  whenever  the  cloud  deck  is  absent.  The  preceding  array  elements 
then  be:ome  the  discriminator  for  determining  the  number  and  type  of 
cloud  decks  present. 

The  broadband  infrared  emlsslvity,  reflectivity,  and  transmissivity 
for  high  clouds  are  functions  of  the  vertical  cloud  liquid  water  content 
following  the  parameterization  equation  developed  by  Liou  and  Wittman 
(1979)  in  the  form 

?  i 

R(W)  =  l  c  W  ,  (3.16) 

i-0 

where  R(W)  denotes  the  flux  reflectivity,  transmissivity,  or  emlsslvity 


3-deck 


LCBOT 


2-deck, 
high  and 
middle  or 
low  clouds 


2-deck, 
middle  and 
low  clouds 


between  clear,  one-,  two-,  and  three-deck  cloudy  atmospheres. 

The  Integral  quantities  LC(1),  LC(3),  and  LC(5)  are  the  cloud  top 
level  indices,  which  are  set  to  zero  if  there  is  no  cloud,  where¬ 
as  the  parameters  LCTOP  and  LCBOT  are  the  high  cloud  top  and 
base  level  indices. 


Table  Predictor  coefficients  ct  for  high  clouds  In  the  Infrared 
region. 

Reflectivity  Transmissivity  Emisslvlty 


0  0.3061 9E-01  0.73597E+00 

1  0.8113^+00  -0. 381 62E+01 

2  -0.1 8995E+01  0.83288E+01 

3  0.24900E+01  -0.91066E+01 

4  -0.15805E+01  0.48765E+01 

5  0.37581E+00  -0. 1 01 41 E+01 


0 . 28568E+00 
0. 35222E+01 
-0.76275E+01 
0.83096E+01 
-0.44405E+01 
0. 92209E+00 


of  the  cirrus  cloud,  c^  the  predictor  coefficients,  and  W  the  vertical 

r>  _p 

cloud  liquid  water  content  in  units  of  10  g  m  .  The  values  of  c^  are 
listed  in  Table  4. 

For  partly  cloudy  conditions,  the  Infrared  cooling  rate  at  each 
model  layer  is  obtained  by  linearly  weighting  the  percentages  of  the 
total  cloud  cover  n  and  clear  portion  (1  -  n)  in  the  form 

-  "  $  +  (1  -  ">  If  •  (3-17) 

*  j 

where  the  superscripts  pc,  c,  and  nc  represent  partly  cloudy,  cloudy,  and 
clear  (no  cloud)  cases,  respectively.  This  equation  is  also  applicable 
to  solar  heating  rate  calculations  under  partly  cloudy  conditions. 

The  radiative  heating  rate  at  level  z  is  related  to  the  divergence 
of  net  fluxes  and  is  given  by 

1  a  1 

at  ‘  -Tcrsr4  <*>  -  F  •  (3-,8) 

P 

where  p  is  the  air  density,  Cp  the  specific  heat  at  constant  pressure, 


and  Az  the  model  layer  thickness.  Equation  (3.18)  is  only  applicable  to 
the  cloudless  region.  To  compute  the  heating  rate  within  clouds,  two 
methods  are  employed.  For  low  and  middle  black  clouds,  an  approximation 
method  has  been  developed.  The  net  flux  at  the  model  level  within  the 
cloud  deck  may  be  evaluated  by  the  difference  between  upward  blackbody 
emission  from  the  lower  layer  and  downward  blackbody  emission  from  the 
upper  layer,  viz., 

F(z)  =  oT4(z-Az/2)  -  0T4(z+Az/2)  .  (3.19) 

Since  T(z-Az/2)  >  T(z+Az/2),  F(z)  is  always  positive.  Near  the  cloud 
top,  since  F(zfc)  =  oT4(zt-Az/2)  -  F+(zt),  and  F'f'(zfc)  «  oT1* (zfc+Az/2) , 
F(zt)  >>  F(zt~Az),  implying  a  strong  cooling  there,  while  near  the  cloud 
base,  F(zb)  >  F(zb+Az),  heating  is  implied.  A  comparison  of  the  result¬ 
ing  Intracloud  cooling  rate  computed  from  this  approximation  with  that 
obtained  from  a  more  detailed  radiative  transfer  model  developed  by  Kinne 
(1987)  will  be  shown  in  a  later  section.  For  high  clouds,  due  to  their 
nonblack  character,  the  above  approximation  cannot  be  applied.  We 
determine  the  intracloud  flux  profile  for  high  clouds  by  exponential 
interpolation  between  cloud  top  and  base  fluxes. 

3.2.3  Solar  Flux  Exchange 

In  our  previous  solar  radiation  program,  calculations  could  only  be 
performed  for  one-  and  two-deck  cloudy  atmospheres.  We  have  expanded  and 
modified  the  program  to  accomodate  the  computation  of  flux  exchanges  in  a 
three-deck  cloudy  atmosphere.  Referring  to  Fig.  4,  we  denote  the  flux 
reflectivity  and  transmissivity  for  the  high  (1),  middle  (2),  and  low  (3) 
clouds  as  (r1  ,t1  ) ,  (r2,t2)  and  (r^.t^),  respectively.  The  associated 


diffuse  quantities  are  defined  by 


Moreover,  the  gaseous  transmissivities, 


u0  dyQ  .  (3.20) 

T0,  T] ,  T2 •  and  T 3  for  respective 


separate  clear  regions  are  defined  by 


where  the  subscripts  0,  1,  2,  and  3  represent  regions  above  the  high 
cloud,  between  high  and  middle  clouds,  between  middle  and  low  clouds,  and 
below  the  low  cloud,  respectively.  The  quantities  AQ,  A1 ,  A2,  and  A^  are 
the  total  absorptivities  of  the  respective  regions.  Moreover,  we  define 
the  multiple  reflection  factor  Sjj  as  the  gain  for  either  the  upward  or 
downward  flux  component  due  to  multiple  reflections  between  the  base  of 
the  1th  cloud  and  top  of  the  jth  cloud  (or  surface).  There  are  three 
multiple  reflection  factors,  which  may  be  expressed  in  the  forms 


where  the  superscript  ♦  signifies  transmission  of  upward  propagating 
radiation,  and  rg  la  the  surface  albedo.  The  over  bars  denote  diffuse 
quantities.  With  these  definitions,  we  then  define  the  generating 
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(3.27) 


-+ 


The  function  G2  Is  also  composed  of  two  components  involving  multiple 


reflections  of  the  direct  downward  flux,  and  transmission  of  the  diffuse 


-t 


upward  flux  through  cloud  2.  Thus, 


% 


a2  r2  S,2  *  S3  *2 


(3.28) 


The  function  G^  is  composed  of  three  components.  These  are  multiple 


reflections  of  direct  and  diffuse  downward  fluxes,  and  upward  transmis- 


-  + 


slon  of  the  diffuse  upward  flux  G^  through  cloud  3.  Thus,  we  write 


G3  =  ^G3  T2  r 3  +  G2  %  *2  T2  r3^  S23  +  G1  ?3  *3  ‘  (3‘29) 


-  + 


The  function  Gjj  is  composed  of  two  components  associated  with  multiple 
reflections  of  direct  and  diffuse  downward  fluxes  at  the  surface,  that 
is, 


5;  ■  («;  n  u  r.  5 


4>  ♦  *•4’  ~  4  —  ^5,4*  ' 

m  ~  rn  »  rn 

3  2  *3  V  ‘s  ^3^ 


(3. 30) 


The  expressions  for  the  downward  diffuse  function  G  may  be  derived  in  a 
similar  manner.  At  the  model  top,  there  is  no  downward  diffuse  radia¬ 
tion,  so  that 


Gj  -  0 


(3-31  ) 


-  i 


Further  downward,  G2  may  be  expressed  as  the  reflected  amount  of  the 
upward  diffuse  flux  in  the  form 


G2  "  G2  T1  rl  ’ 


(3-32) 
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a 


i 


i 
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a 
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Subsequently,  5^  and  may  be  wnltten  as  the  sum  of  the  reflected  amount 
of  both  diffuse  upward  and  diffuse  downward  fluxes  In  the  forms 


s;  •  e;  n  *  %  n  ^ 


(3.33) 


5I-5IT3?3+  *3^*3  • 


(3.31*) 


4-  4 

The  G  functions  relevant  to  direct  downward  transmission,  i.e.,  G-j  ,  Gp, 
■4  4 

G^ ,  and  G^ ,  Involve  known  quantities  of  cloud  and  gaseous  absorption. 

4  4 

They  can  be  successively  computed  explicitly  In  the  order  of  G1  -»  C2  ♦ 

4  4 

G3  -*■  Gi| .  The  rest  of  the  G  functions  are  interelated.  By  substituting 
Eq,  (3-28)  into  Eq.  (3-32),  we  obtain  G2  as  a  function  of  G^: 


°2  '  1°)  <  r2  St2  *  °3  T2  €2>  F1  '  (3’35) 

Moreover,  substituting  Eq.  (3.35)  into  Eq.  (3.30),  leads  to  an  expression 
relating  G^  as  a  function  of  G^  In  the  form 


[«; *  c;  -a  z  v  F< 


?;  t3  t\\ 


r 

s 


Finally,  we  substitute  both  Eqs.  (3.35)  and  (3.36)  into  Eq. 
solve  for  G^  to  obtain 


(3.36) 
(3.29),  ard 


0*  "  (y1  +  v2  +  y3  +  yJ  '  ^  ‘  X1  "  x2 ^ 


(3.37) 


where 


g3  T2  r3  S23 


(3.38) 
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(3.39) 


G2  r2  S12  F1  X  *2  %  r3  S23  * 


G4  T3  rs  S34  T3  S  ’ 


(3.40) 


G2  T1  r2  S12  T\  rl  ^1  l2  T2  t3  T3  rs  S34  "^3  ^  ’  (3*4l) 


T2  fc2  T1  rl  T1  l2  T2  r3  S23  ’ 


(3.42) 


X2  =  ^2  r2  +  T2  t2  T\  P1  ri  fc2^  T2  t3  T3  rs  S34  T3  ^3  * 

(3.43) 

The  remaining  G  functions  may  be  computed  according  to  Eqs.  ( 3 . 27 )— (3.34) 
in  the  order  V2  +  g|  -►  G2  -»  G3  ■»  GjJ  -  G4. 

Using  the  preceding  G  functions,  we  can  compute  the  flux  components 
for  each  clear  region  separated  by  cloud  decks.  Above  the  high  cloud  we 


F+(z)  =  u0  SQ  Q*  [l  -  A+(zt-z)]  , 


(3.44) 


F*(z)  «  p0  SQ  Gj  [1  -  Af(z-2t1)J  . 


(3-45) 


where  SQ  denotes  the  solar  constant,  and  the  cosine  of  the  solar 
zenith  angle.  Between  the  high  and  middle  clouds,  we  find 


F*(z)  =  uQ  SQ  lG*  [l  -  A*(zb1-z)]  ♦  G*  [l  -  A+(zb1-z)](  ,  (3.46) 


F  (z) 


•  "0  S0  C2  U  •  A,(z-zt2))  • 


(3.47) 


Between  the  middle  and  low  clouds,  we  have 
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F+(z)  -  mq  SQ  (g*  [l  -  A*(zb2-z)]  +  G*  [ 1  -  A+(zb2-z)]}  ,  (3-^8) 


F+ (z)  -  uQ  SQ  G*  [1  -  A*(z-zt3)]  .  (3.49) 


And  below  the  low  cloud,  we  find 


F* (z)  «  uQ  S0  {G*  [1  -  A+(zb3-z)]  +  G*  [1  -  A+(zb3-z)]}  ,  (3.50) 

l 

l 

F+(z)  -  u0  SQ  Gj  [l  -  A + ( z )  ]  .  (3.51) 


The  broadband  solar  reflection  and  transmission  values  for  various 
cloud  decks  are  obtained  through  the  following  parameterization  equation 
based  on  detailed  radiative  transfer  calculations: 


r.t(w0,W) 


3  3 

l  l 


m*0  n«0 


b 


mn 


(3.52) 


where  bmn  are  empirical  coefficients  derived  from  the  parameterization 
principle  proposed  by  Llou  and  Wlttman  (1979),  and  the  vertical  liquid 
water  content  W  is  In  units  of  102  g  m-2.  The  bmn  for  Cu  (low  cloud), 

As  (middle  cloud),  and  Ci  (high  cloud)  are  given  in  Tables  5,  6,  and  7, 
respectively.  Computations  of  solar  heating  rates  in  clear,  overcast, 
and  partly  cloudy  conditions  follow  Eqs.  (3.17)  and  (3.18).  However,  to 
compute  intracloud  heating  rates,  we  use  an  exponential  interpolation  of 
fluxes  at  the  cloud  top  and  base.  Results  will  be  presented  in  the  next 
section . 


32 


■'V'. ' *■ .  A  AO  A  AAA  .\A  AAA  A.  ■ A.  AA.\  .' 

iwt.  ifc.  w1-  j ^ -  A  Wl.’jL.  *-*  d*...  *-*  S  ^  W  '-ik'j.  V 


BUTWH  VIM  W  UIU  U V.1 ".ur LWKV  V  V  V  vv 


w.v.v.v  v.v  v  v  v  v  v.vv.v.v.v.'^V.’/.VWW' 


L 


Table  5.  Approximating  predictor  coefficients  bjj  for  Cu  (low  cloud). 


j 

1 

0 

1 

2 

3 

Reflection 

0 

0 . 51 771 E+00 

0.18277E+00 

-0 . 35851 E  —0 1 

0.23i478E-02 

1 

0.72387E+00 

-0.2701 9E+00 

0.51 303E-01 

-0. 3^^28E-02 

2 

-0.21 H99E+01 

0.75225E+00 

-0.14588E+00 

0.97306E-02 

3 

0. 1 21 57E+01 

-0.1J0120E+00 

0.781  M5E-01 

-0.52224E-02 

Transmission 

0 

0.i)3642E+00 

-0.22727E+00 

0.U4702E-01 

-0.291 1 7E-02 

1 

-0.6631 9E+00 

0.35^1 5E+00 

-0.70708E-01 

0.146403E-02 

2 

0.19014E+01 

-0.10256E+01 

0.2057HE+00 

-0. 1  353  ^E— 01 

3 

-0 . 1 01 57E+01 

0.5i1533E+00 

-0.1 091 2E+00 

0. 71 667E-02 

Absorption 

0 

0.66103E-01 

0. 3^702E-01 

-0.69662E-02 

0.W253E-03 

1 

-0.79148HE-01 

-0. 7361 4E-01 

0. 1  7572E-01 

-0.10799E-02 

2 

0.29293E+00 

0.25047E+00 

-0.55398E-01 

0.351 39E-02 

3 

-0.22077E+00 

-0. 1 3272E+00 

0.28740E-01 

-0.17990E-02 

Table  6.  Approximating  predictor  coefficients  bjj  for  As  (middle 
cloud ) . 
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v«r 


9 


1 


a 


i 


1 

0 

1 

2 

3 

r 

*5 

y 

& 

Reflection 

0 

0.6U1457E+00 

0. 28785E+00 

-0. 1 2069E+00 

0.1 6673E-01 

1 

0.U14  388E+00 

-0.il51i96E+00 

0.1 8695E+00 

-0.26672E-01 

$ 

2 

-0. 1  14875E+01 

0. 1  31  *41  E+01 

-0.54825E+00 

0.771 21 E-01 

3 

0.85585E+00 

-0.71 3^6E+00  0.29823E+00 

Transmission 

-0.141  922E-01 

% 

A 

0 

0. 30735E+00 

-0. 35373E+00 

0.1  14797E+00 

-0.20208E-01 

1 

-0.11H042E+00 

0  50188E+00 

-0.20975E+00 

0.28630E-01 

yj 

2 

0.1 3537E+01 

-0 . 1 5738E+01 

0.662H8E+00 

-0.90739E-01 

3 

-0. 73396E+00 

0. 85256E+00 

Absorption 

-0. 35858E+00 

0.149077E-01 

1 

| 

0 

0 . 705914E-01 

0.no901E-01 

-0.16982E-01 

0.21  3145E-02 

1 

-0.286M0E-01 

-0 . 1 891 7E-01 

0.111 67E-01 

-0.37025E-03 

* 

2 

0. 1 9500E+00 

0. 1 9298E+00 

-0.86520E-01 

0.98269E-02 

£ 

3 

-0.151 33E+00 

-0.10728E+00 

0.1471 143E-01 

-0.535214E-02 

c 

a 


Table  7.  Approximating  predictor  coefficients  bj_j  for  Cl  (high  cloud). 


j 

1 

0 

1 

2 

3 

0 

0.25264E+00 

Reflection 

0.45462E+00 

-0.29492E+00 

0.82536E-01 

1 

-0.29273E+00 

0. 22029E+01 

-0.24834E+01 

0.82674E+00 

2 

-0.2U531E+00 

-0.63610E+01 

0.75864E+01 

-0.25896E+01 

3 

0. 35097E+00 

0.38933E+01 

-0.49188E+01 

0.171 87E+01 

0 

0.78722E+00 

Transmission 

-0.66348E+00  0.44495E+00 

-0 . 1 2500E+00 

1 

0.73255E+00 

-0.43772E+01 

0.47094E+01 

-0 . 1 5549E+01 

2 

-0.68508E+00 

0.11764E+02 

-0.13716E+02 

0. 46796E+01 

3 

0. 1 5453E+00 

-0.70865E+01 

0.87472E+01 

-0.30566E+01 

0 

0.49388E-01 

Absorption 
0. 1 371 5E+00 

-0. 1 01 22E+00 

0.28402E-01 

1 

-0. 1 9622E+00 

0. 1 1882E+01 

-0. 12267E+01 

0. 40299E+00 

2 

0.241 18E+00 

-0.231 49E+01 

0.28776E+01 

-0 . 1 01 68E+01 

3 

-0.94420E-01 

0.11 305E+01 

“0. 1  5772E+01 

0.58577E+00 

Section  M 


INTERCOMPARISON  OF  RADIATIVE  HEATING 
RATE  PROFILES 

In  this  section,  we  present  heating  rate  results  for  clear  and 
various  cloudy  atmospheres. 

4 . 1  Intercomparison  of  Clear-Sky  IR  Heating  Rates 

Using  the  IR  radiative  transfer  parameterization  described  in 
Section  2,  we  have  carried  out  a  number  of  numerical  comparisons  and 
verifications  with  more  detailed  band  model  and/or  line-by-line 
calculations  (Liou  and  Ou,  1981;  Ou  and  Liou,  1983).  Specifically,  we 
have  computed  cooling  rates  due  to  H20  and  0^  emission,  in  a  tropical 
atmosphere  and  compared  them  with  results  from  a  band  model  developed  by 
Roewe  and  Liou  (1978),  and  the  line-by-line  calculations  reported  by  Chou 
and  Arklrg  (1980).  Differences  between  these  sets  of  calculations  are 
within  about  0.5°K  day-1.  In  addition,  cooling  rates  due  to  C02  emission 
in  both  standard  and  subarctic  atmospheres  have  been  compared  with 
results  from  line-by-line  computations  presented  by  Fels  and  Schwarzkopf 
(1981).  The  difference  near  maximum  cooling  ('50  km)  is  no  more  thant  5% 
(~0.1i°K).  The  total  IR  cooling  rate  profile  for  a  tropical  atmosphere 
also  agrees  closely  with  that  computed  from  a  band  model  (Roewe  and  Liou, 
1978)  and  from  observations  reported  by  Cox  (1969). 

In  addition  to  the  aforementioned  comparisons,  we  perform  computa¬ 
tions  of  the  clear-sky  IR  radiative  cooling  rate  profiles  for  five 


different  model  atmospheres  provided  by  McClatchey  et  al .  (1971),  and 
compare  them  with  results  produced  by  the  GFDL/NMC  radiation  parameter¬ 
ization  method.  The  current  model  computations  are  based  on  the  AFGL  18 
layer  vertical  grid  structure,  which  has  been  shown  in  Fig.  5.  The  AFGL 
18-layer  structure  is  identical  to  the  NMC  18-layer  structure.  The  AFGL 
1 8-layer  model  has  four  layers  in  the  stratosphere  (o  <  0.2)  and  four 
layers  in  the  boundary  layer  region  (o  >  0.875).  The  cooling  rate 
profiles  are  shown  in  Fig.  6a-e.  In  general,  the  IR  cooling  rate  pro¬ 
files  computed  from  the  present  model  agree  well  with  the  GFDL/NMC 
results.  The  differences  between  these  two  models  are  generally  within 
about  0.5°K/day,  except  in  the  lower  atmosphere.  The  largest  difference 
is  found  in  the  tropical  lower  atmosphere,  probably  due  to  the  treatment 
of  water  vapor  continuum  absorption. 

14 .2  Adding  Scheme  for  Radiative  Transfer  in  Cloudy  Atmospheres 

A  detailed  narrow-band  radiative  transfer  program  was  recently 
developed  by  Klnne  (1987).  This  program  employs  the  matrix-operator 
method  (^lass  et  al . ,  1973),  which  is  essentially  the  same  as  the  adding 
method  in  terms  of  numerical  computations.  The  solar  and  IR  spectra  are 
divided  into  a  number  of  bands.  In  each  band,  the  transmittance  of 
atmospheric  gases  is  approximated  by  a  set  of  exponential  functions. 

Application  of  the  adding  principle  to  an  inhomogeneous  atmosphere 
is  carried  out  by  dividing  the  atmosphere  into  a  number  of  homogeneous 
layers.  Two  steps  are  followed.  First,  the  reflection  R,  transmission 
,  and  emission  J  of  a  layer  are  determined.  A  very  thin  sublayer  with 
an  optical  depth  t  <  10~^  is  used  initially.  For  these  sublayers, 


Fig.  6  Intercomparison  of  model-generated  clear-sky  IR  cooling  rate  profiles  with  those  from  the 

GFDL/NMC  model  for  (a)  tropical,  (b)  midlatitude  summer,  (c)  midlatitude  winter,  (d)  subarctic 
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matrices  R  and  T  and  the  vector  J  are  determined  In  a  manner  such  that 
the  effects  of  direct  solar  radiation  and  gaseous  absorption  are 
Included.  Then  two  layers  with  the  same  optical  properties  are  combined, 
until  the  desired  optical  thickness  for  a  prescribed  layer  is  reached. 
This  Is  the  doubling  procedure,  which  is  performed  for  all  prescribed 
homogeneous  layers. 

There  is  now  a  set  of  R,  T,  and  J  for  all  homogeneous  layers.  Two 
layers  are  subsequently  added.  Beginning  from  the  top  of  the  atmosphere, 
using  the  optical  properties  of  the  uppermost  layer,  the  values  of  pro¬ 
gressively  lower  layers  are  added  one  at  a  time.  This  adding  procedure 
Is  applied  for  all  homogeneous  layers.  Finally,  a  surface  layer  Is  added 
to  account  for  the  emission  and  reflection  effects  of  the  surface. 

The  resulting  Intensities  I  for  various  directions  at  each  layer 

t  4- 

boundary  are  summed  to  obtain  the  hemispheric  fluxes  F  and  F  .  From  the 
divergence  of  the  net  flux  FN  at  the  boundaries  of  a  given  layer,  the 
layer  heating  and  cooling  rates  are  determined.  The  cooling  rate  results 
to  be  shown  In  this  section  were  obtained  based  on  a  spectral  resolution 
of  8  solar  and  12  IR  bands  approximated  by  exponential  functions  of  a 
total  of  120  terms,  and  on  five  equal  intervals  of  the  zenith  angle.  The 
Inhomogeneous  atmosphere  Is  divided  Into  18  layers. 

With  the  exception  of  the  ozone  data  at  solar  wavelengths,  which  are 
taken  from  measurements  by  Vlgroux  (1953),  all  data  for  gaseous  absorp¬ 
tion  are  based  on  spectral  transmittance  data  computed  from  line-by-llne 
data  compiled  by  McClatchey  et  al .  (1971)  and  updated  by  Rothman  et  al. 
(1983).  For  water  clouds,  the  spectral  distribution  of  refractive 
Indices  Is  taken  from  values  provided  by  Downing  and  Williams  (1975),  and 


the  single-scattering  albedo  and  asymmetry  factor  were  computed  using  Mie 
scattering  theory.  For  ice  clouds,  the  refractive  indices  for  ice  were 
taken  from  values  tabulated  by  Schaaf  and  Williams  (1973).  The  single¬ 
scattering  albedo  and  asymmetry  factor  were  approximated  by  those  of 
spherical  water  droplets  with  the  same  surface  area,  but  with  adjustments 
to  account  for  the  hexagonal  structure  of  ice  particles.  The  particle 
size  distributions  used  for  single-scattering  calculations  for  high, 
middle,  and  low  clouds  follow,  respectively,  those  for  cumulus,  alto- 
stratus,  and  cirrostratus  (Liou  and  Wittman,  1979). 

ij.3  Intercomparison  of  Cloudy-Sky  Radiative  Heating  Rates 

Utilizing  the  new  scheme  for  computing  the  radiative  heating  rate  in 
cloudy  atmospheres  as  described  in  Section  3,  we  carry  out  intercom¬ 
parisons  between  results  from  the  present  model  and  those  from  the 
doubling/adding  scheme.  The  present  computation  is  based  on  the  AFGL  18- 
layer  vertical  grid  structure,  as  shown  in  Fig.  5.  Figure  5  also  shows 
the  vertical  domain  for  each  of  the  three  cloud  decks.  Clouds  are  assumed 
to  fill  respective  vertical  domains.  The  atmospheric  profile  used  for 
flux  calculations  is  the  1962  U.S.  Standard  Atmosphere.  Computations 
were  made  for  seven  cloud  configurations  including  three  one-deck,  three 
two-deck,  and  one  three-deck  clouds.  In  the  infrared  region,  the  high 
cloud  radiative  properties  were  prescribed  based  on  a  liquid  water  con¬ 
tent  of  0.006  g/m^.  in  the  solar  region,  the  cloud  radiative  properties 
were  calculated  using  the  parameterlzations  of  Liou  and  Wittman  (1979). 
The  cosine  of  the  solar  zenith  angle  is  chosen  to  be  0.8  and  the  surface 


albedo  is  0.2. 


Figure  7a-h  shows  the  results  of  IR  radiation  intercoinpar isons .  The 


plots  include  IR  radiative  cooling  rates  obtained  from  the  current 
parameterization  method,  as  well  as  from  the  doubling/adding  model  using 
the  AFGL  18-layer  structure.  In  general,  the  cooling  rate  profiles  agree 
with  each  other  rather  well  in  the  clear  atmosphere  case  with  deviations 
less  thar  0.3°C/day  (Fig.  7a).  They  are  also  quite  close  to  each  other 
for  clear  regions  separated  by  clouds  as  shown  in  Fig.  7b-h.  Larger 
differences  are  noted  between  intracloud  cooling  rate  values.  For  high 
clouds,  agreement  between  the  current  parameteri zation  and  the  doubling/ 
adding  scheme  is  better  near  the  cloud  top  and  base  than  in  mid-cloud 
regions,  where  differences  are  on  the  order  of  -1°K/day.  The  present 
parameteri zation  reproduces  quite  realistically  cloud  top  cooling  and 
base  warming.  This  appears  to  demonstrate  the  validity  of  the  exponen¬ 
tial  interpolation  for  computations  of  Intracloud  fluxes.  For  middle  and 
low  clouds,  excellent  agreement  is  shown  between  the  two  sets  of  results, 
suggesting  the  propriety  of  the  blackbody  assumption.  The  cloud  top 
cooling  for  middle  and  low  clouds  is  on  the  order  of  -10  and  -12°K/day, 
respectively.  In  the  cases  with  two-deck  clouds,  cooling  is  reduced  due 
to  downward  emission  by  the  higher  cloud  base.  The  cloud  base  heating 
for  middle  and  low  clouds  is  on  the  order  of  3.5  and  1.7°K/day,  respec¬ 
tively.  In  mid-cloud  regions,  there  is  no  heating  as  a  result  of 
constant  net  fluxes,  which  are  close  to  zero.  Zero  net  flux  indicates 
that  the  upward  and  downward  fluxes  at  that  level  are  equal.  In  the  case 
wher  both  middle  and  low  clouds  are  present  (Fig.  7g,h),  cooling  near  the 
boundary  of  these  two  clouds  is  also  zero,  sirce  the  two  cloud  decks  are 


contiguous . 
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Figure  8a-h  shows  the  results  of  solar  radiatior  irtercompar isors . 

In  general,  heating  rate  profiles  for  the  clear  atmosphere  agree  well 
with  each  other  in  the  troposphere,  with  differences  of  no  more  than 
0.5°K/day.  In  the  stratosphere,  the  present  parameter  1 zat ion  produces 
slightly  larger  heating  rates  than  the  doubl ing/addirg  scheme.  Good 
agreements  are  also  evident  for  heating  rate  pnoflles  in  clear  regiors 
separated  by  clouds,  as  shown  in  Fig.  8b-h.  Withir  each  cloud  deck, 
larger  differences  are  shown.  These  differences  may  be  caused  by  the 
use  of  broadband  cloud  radiative  properties  in  the  parameterization.  Ir 
general,  the  present  parameterization  cornectly  produces  the  pattern  of 
cloud  top  heatirg.  For  high  clouds,  the  parameterization  produces  less 
heating  than  the  doubling/adding  scheme,  but  the  difference  is  no  more 
thar  1°K/day.  For  middle  and  low  clouds,  the  difference  between  the 
current  parameterization  and  the  doubling/addirg  scheme  is  within 
2°K/day.  In  the  case  of  one-deck  clouds  (Fig.  8b-d),  heating  peaks  of 
about  8  and  9.5°K/day  at  the  cloud  top  are  generated  by  middle  and  low 
clouds,  respectively.  In  the  case  of  two-deck  clouds  (Fig.  8e-g),  if  the 
uppen  deck  is  the  high  cloud,  the  heating  rate  profiles  for  middle  or  low 
clouds  do  not  change  signif Icartly .  However,  for  a  combination  of  middle 
and  low  clouds,  heatirg  rates  below  the  optically  thick  middle  cloud 
become  very  small.  This  is  also  evident  for  the  three-deck  cloud  case. 

Ir  summary,  the  presert  broadband  parameterization  programs  for  the 
transfer  of  IR  and  solar  fluxes  produce  reliable  heating/cooling  profiles 
ir  various  cloudy  conditions.  The  maximum  deviation  of  about  2°K/day 
occurs  for  the  solar  heatirg  rate  at  the  low  cloud  top.  Other  deviations 
produced  by  this  efficiert  ard  economic  parameter! zat ior  scheme  are  all 


withir  about  1°K/day. 


Intercompar ison  of  model-generated  solar  heating  rate  profiles  with  those  from  the  doubling/ 
ldding  scheme  based  on  the  1962  U.S.  Standard  Atmosphere  for  different  cases:  (a)  clear, 

(b)  high  cloud  only,  (c)  middle  cloud  only,  (d)  low  cloud  only,  (e)  high  and  middle  clouds, 
(f)  high  and  low  clouds,  (g)  middle  and  low  clouds,  and  (h)  high,  middle,  and  low  clouds. 


Section  5 


RESULTS  OF  CLOUD  AND  RADIATION  BUDGET  FIELDS 

5.1  Cloud  Prediction 

Using  the  AFGL  12-layer  GCM  and  the  Initial  field  of  12Z,  17  June 
1979,  which  consists  of  temperature  and  specific  humidity  data,  clouds 
are  formed  according  to  the  method  described  In  subsection  3.1.1.  The 
multiple-layer  cloud  field  is  then  strapped  into  at  most  three  decks  of 
clouds  using  the  scheme  described  In  subsection  3.1.2.  Subsequently,  the 
total  cloud  cover  field  is  generated  using  the  maximum  overlap  assump¬ 
tion.  The  zonally-averaged  latitudinal  distribution  of  high,  middle, 
low,  and  total  cloud  cover  for  one  time  step  Is  shown  in  Fig.  9a-d.  The 
high  cloud  cover  Is  mostly  -10$,  except  near  the  South  Pole,  where  It 
reaches  more  than  40$.  The  middle  cloud  cover  reaches  ^0$  near  10°N  and 
80$  near  80°S,  but  decreases  to  less  than  10$  In  the  subtopics.  The  low 
cloud  distribution  Is  similar  to  that  for  middle  clouds  except  near  the 
North  Pole  region,  where  more  low  clouds  are  present.  The  total  cloud 
cover  distribution  Is  slightly  lower  than  the  summer  cloud  climatology  of 
London  (1957). 

5.2  Radiation  Budget 

The  global  distributions  of  outgoing  IR  and  net  downward  solar 
fluxes  at  the  top  of  the  atmosphere  (T0A)  and  surface  are  computed  from 
the  AFGL  12-layer  GCM.  The  zonally  averaged  latitudinal  distributions  of 
these  quantities  are  shown  In  Fig.  lOa-d.  The  outgoing  IR  flux  at  the 


CLOUD  COVER  (%) 


i'i  j't  t>a'«-i1«a  ».>  '»ki  ViViV.* 


(b)  MIDDLE  CLOUD 


(d)  TOTAL  — °~  ~  LO NDO 

-Os  -  MODEL 


>o  ' 
""</  / 


LATITUDE  (deg) 


Q  Zonal ly  aviraged  latitudinal  clou' 

the  Initial  temperature  and  humid! 
(a)  high  cloud,  (b)  middle  cloud, 


TOA  (Fig.  10a)  is  lowest  near  the  South  Pole  (-110  W/m2),  and  increases 
northward  to  about  260  W/m2  in  the  southern  subtropics.  It  decreases  to 

240  W/m2  and  Increases  to  270  W/m2  near  30°N.  An  IR  flux  of  200  W/m2  is 

found  near  the  North  Pole.  Low  values  in  the  Southern  Hemisphere  are 
caused  by  the  low  temperatures  of  the  southern  winter.  The  minimum  in 
the  tropics  is  due  to  the  combined  effects  of  high  humidity  and  cloudi¬ 
ness.  The  net  solar  flux  at  the  TOA  (Fig.  10b)  is  0  near  the  South  Pole, 
increases  to  280  W/m  near  5°S  before  a  slight  decrease  to  230  W/m  .  It 

increases  to  a  maximum  of  -400  W/m2  near  35°N  and  drops  to'  230  W/m2  near 

the  North  Pole.  The  maximum  in  the  northern  subtropics  is  due  to  the 
combined  effects  of  high  solar  zenith  angles  in  the  polar  region  and 
large  fractional  sunlit  areas  in  this  latitudinal  belt.  The  peak  value 
(-400  W/m2)  is  larger  than  the  satellite-derived  maximum  (-340  W/m2)  for 
the  summer  season  (Stephens  et  al . ,  1981).  This  is  expected  since  the 
present  calculation  is  for  17  June,  which  is  only  five  days  away  from  the 
summer  solstice.  The  distribution  of  surface  net  IR  fluxes  (Fig.  10c) 
shows  that  the  average  surface  IR  flux  varies  between  20  and  90  W/m2.  A 
large  dip  is  present  in  the  equatorial  area,  due  to  the  combined  effect 
of  clouds  and  humidity.  Finally,  the  pattern  of  net  downward  solar 
fluxes  at  the  surface  is  similar  to  that  of  net  solar  fluxes  at  the  TOA, 
except  the  values  are  smaller  by  about  0-100  W/m2.  This  difference  is 
attributed  to  atmospheric  absorption,  which  depends  on  the  solar  zenith 
angle  and  humidity  distribution.  The  highest  absorption  occurs  near  the 
equatorial  region  due  to  higher  humidity. 

5. 3  Radiative  Cooling  Rates  in  the  Atmosphere 

Net  radiative  heating  rates  were  calculated  using  the  current 
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radiation  and  cloud  model  with  a  12-layer  structure.  The  zonally 
averaged  latitude-sigma  distributions  of  the  net  radiative  heating  rate 
for  clear  and  model-generated  cloudy  atmospheres,  as  well  as  the 
difference  between  the  two  are  presented  in  Fig.  Ila-c.  For  the  clear 
case  (Fig.  11a),  the  net  cooling  in  the  troposphere  increases  from  the 
tropopause  toward  the  surface.  It  ranges  between  0  and  -4°K/day.  The 
large  cooling  near  the  surface  layers  in  the  equatorial  region  is  due  to 
the  contribution  of  water  vapor  continuum  absorption.  Somewhat  larger 
cooling  is  present  near  the  South  Pole  due  to  the  lack  of  solar  insola¬ 
tion  in  the  antarctic.  Near  the  North  Pole,  however,  there  is  a  slight 
warming  due  to  absorption  of  solar  radiation  in  the  long  arctic  day.  The 
northern  stratosphere  also  shows  a  slight  warming  due  to  ozone  absorp¬ 
tion. 

Compared  with  the  clear  case,  the  net  cooling  produced  by  cloudy 
atmospheres  is  smaller  in  the  tropical  troposphere  and  the  lower  tropo¬ 
sphere  of  the  Southern  Hemisphere.  The  effect  of  clouds  on  the  net  cool¬ 
ing  can  be  seen  in  Fig.  11c,  where  the  latitude-sigma  distribution  of 
net  cooling  differences  is  shown.  Below  o  =  0.84,  clouds  have  a  warming 
effect.  Between  o  *  0.64  and  0.84,  several  areas  of  cooling  are  present, 
but  cooling  is  not  present  between  the  equator  and  60°N,  presumably  due 
to  absorption  of  strong  solar  solar  Insolation  by  clouds.  Above  o  = 

0.51,  clouds  generate  cooling  everywhere  in  the  uppen  troposphere.  Since 
the  model  middle  and  high  cloud  tops  are  located  at  o  *=  0.45  and  0.2, 
respectively  (Fig.  5),  the  net  cooling  there  is  directly  associated  with 
strong  cloud  top  cooling.  In  addition,  the  IR  emission  of  high  clouds 
outweighs  the  absorption  due  to  solar  radiation.  As  a  result,  the 
presence  of  high  clouds  produces  cooling  at  the  cloud  top. 


(c)  DIFFERENCE  (°K/DAY) 
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Section  6 

SUMMARY 

In  this  report,  we  have  described  the  IR  and  solar  radiation  param¬ 
eterization  programs  for  clear  and  generalized  three-deck  cloudy  atmo¬ 
spheres  for  incorporation  in  the  AFGL  GCM.  For  the  IR  part,  high  clouds 
are  considered  to  be  nonblack,  but  middle  and  low  clouds  are  assumed  to 
be  blackbodles.  The  flux  profiles  in  clear  regions  separated  by  clouds 
and/or  the  surface  are  computed  based  on  the  same  principle  used  in 
computing  clear  atmosphere  fluxes.  Fluxes  within  the  clouds  are  deter¬ 
mined  by  interpolation  for  nonblack  clouds  and  by  the  blackbody  emission 
for  black  clouds.  For  the  solar  part,  there  is  no  distinction  between 
black  and  nonblack  clouds.  Instead,  cloud  radiative  properties  are 
obtained  based  on  the  prescribed  vertical  cloud  liquid  water  content. 
Fluxes  between  cloud  decks  are  obtained  by  a  ray  tracing  method,  as 
described  in  Section  3. 

Prediction  of  cloud  cover  and  cloud  configurations  is  carried  out 
in  three  steps.  First,  the  cloud  cover  is  computed  for  each  of  the  model 
layers  designated  for  cloud  formation.  Then,  by  strapping  two  or  more 
contiguous  layers  into  cloud  decks,  the  cloud  cover  for  each  deck  is 
computed  based  on  the  averaging  scheme.  Finally,  the  total  cloud  cover 
is  computed  by  the  maximum  overlap  assumption.  There  are  seven  cloud 
configurations,  including  three  one-deck  cases,  three  two-deck  cases,  and 
one  one-deck  case. 

Using  five  atmospheric  profiles,  we  show  that  IR  clear-sky  cooling 
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rates  computed  from  the  present  model  agree  with  those  obtained  from  the 
GFDL/NMC  model  within  about  0.5°K/day.  One  exception  is  near  the 
surface,  where  a  difference  of  about  2°K/day  is  found  for  the  moist 
tropical  atmosphere.  This  difference  is  probably  due  to  different  treat¬ 
ments  of  the  water  vapor  continuum  absorption.  We  have  also  compared 
heating  rates  for  various  cloudy  conditions  computed  from  the  present 
parameterization  program  with  those  from  a  more  exact  radiative  transfer 
method.  This  method  involves  the  use  of  adding/doubling  principles  and 
numerous  spectral  intervals  covering  solar  and  IR  spectra.  We  show  that 
the  present  broadband  parameterization  program,  which  is  computationally 
efficient,  produces  reasonable  solar  heating  and  IR  cooling  rate  pro¬ 
files.  We  have  performed  comparisons  for  atmospheres  containing  one, 
two,  and  three  cloud  decks.  The  largest  deviation  occurs  when  nonblack 
and  semitransparent  high  clouds  are  present. 

The  present  radiation  and  cloud  parameterization  program  has  been 
successfully  Incorporated  in  the  AFGL  12-layer  GCM.  Simulations  have 
been  performed  for  cloud  and  radiation  budget  distributions  using  the 
initial  humidity  and  temperature  fields  corresponding  to  12Z,  17  June 
1979.  Zonally  averaged  global  distributions  of  cloud  cover  and  radiation 
budgets  at  the  TOA  and  surface  are  displayed  and  examined.  Cloud  cover 
and  IR  and  solar  fluxes  at  the  TOA  are  in  general  agreement  with  avail¬ 
able  cloud  and  radiation  budget  climatologies.  In  addition,  we  also 
present  the  zonally  averaged  atmospheric  cooling  rate  profiles  with  and 
without  the  presence  of  cloud  fields  using  the  AFGL  12-layer  GCM.  The 
importance  of  clouds  in.  the  generation  of  atmospheric  cooling  is  pointed 
out . 
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Appendix  A 

SIMPLIFICATION  OF  RADIATIVE  TRANSFER 
IN  HIGH  CLOUDS 


The  upward  and  downward  fluxes  below  the  high  cloud  in  our  original 
program  (Liou  and  Ou,  1981;  Liou  et  al.,  198M)  are  given  by 


F*(z)  -  aT4  [l  -  ef(z,Tj]  +  [  cT4(z’)  K(z-z* )  dz'  , 

s  3  Jq 


( A .  1  ) 


F*(z)  -  rC  (oT4  [l  -  eI(2zb-z,Ta]]  ♦  j  oT^z')  K(2zb~z-z')  dz' 


+  tc  lj  1  oT4(z')  K(z'-zt+zb-z)  dz'}  +  e°  oTb  [l  -  e1  (zb-z.Tb) ] 


c  ^ 


oT  (z')  K(z'-z)  dz'  , 


(A. 2) 


where  the  subscripts  t  and  b  denote  the  high  cloud  top  and  base,  respec¬ 
tively.  For  the  region  above  the  high  cloud,  we  have 


r  Z 

F+(z)  -  rC  oT4(z')  K(z'+z-2z  )  dz'  ♦  tC  |oT4  [l  -  ef( z-zt+zb,Tg] ] 


f  z 

♦  j  b  oT4(z’)  K(z-zt+zb-z')  dz'}  *  eC  oT4  [l  -  ef(z-zt,Tt)] 


+|  oT  { z ' )  K(z-z' )  dz'  , 
z^ 


(A. 3) 


F+(z) 


oT  (z' )  K(z'-z)  dz'  . 


(A.H) 


6** 


In  Eqs.  (A.I)-(A.^)  all  the  notations  have  been  defined  In  Section  2.  As 


shown  In  the  equations,  computations  of  the  broadband  flux  emlsslvity 
become  quite  Involved.  To  economize  the  computational  effort,  we  first 
compute  flux  components  above  and  below  a  black  cloud  at  the  same  posi¬ 
tion  as  the  high  cloud.  Thus,  below  and  above  the  black  cloud,  we  have 


FJ(z) 


Fb(2) 


F*(z) 

0 


1 2 

oT^  [l  -  ef(z,Ts)]  +  |  oAz')  K(z-z')  dz'  , 
ot£  [l  -  ef(zb-z,Tb)]  +  f  oAz')  K(z'-z)  dz’  , 

*  *7 


(A. 5) 


(A. 6) 


4  f*  fz  li 

oT  [l  -  e  (z-z  ,T  )]  +  oT  (z»)  K(z-z’)  dz'  ,  (A. 7) 


f;(z) 


oT  (z')  K(z'-z)  dz'  , 


(A. 8) 


where  the  subscript  b  represents  the  black  cloud.  For  below  the  high 
cloud,  by  substituting  Eqs.  (A. 6)  and  (A. 7)  into  Eqs.  (A. 2)  and  (A. 3), 
respectively,  we  obtain, 


F+(z)  -  F*(z)  +  rC  F* (z)  +  tC  F*(z) 

-  (l  -  eC)  oTb  [l  -  ef(zb-z,Tb)j  ,  (A. 9) 

and  for  above  the  high  cloud, 

F+(z)  =  F*(z)  +  rC  F*(z)  +  t°  F*(z) 
b  r  t 

-  (1  -  eC)  otJ  [l  -  ef(z-zt,Tt)]  ,  (A. 10) 

where  rc  Fp(z),  tc  F^(z),  rc  F,t(z),  and  tc  F^(z)  are  the  reflected  and 


transmitted  flux  components  below  and  above  the  cloud,  respectively.  We 


now  approximate  these  four  quantities  as  follows: 

rC  F*(z)  =  rC  F+(zb)  [l  -  ef(zb-z,Tb)]  ,  (A. 11) 
tC  F*(z)  ■=  t°  F*(zt]  [l  -  ef(zb-z,Tb)]  ,  (A. 12) 
rC  F*(z)  -  rC  F+(zfc)  [l  -  ef(z-zt ,Tt ) ]  ,  (A. 13) 
t°  F*(z)  =  tC  F+(zb]  [l  -  ef(z-zt,Tt)]  ,  (A. 14) 


4,  f 

where  F  (zt)  and  F  (zb)  denote  radiative  fluxes  entering  the  cloud  at  the 
cloud  top  and  bottom,  respectively.  Substituting  Eqs.  (A. 11 ) — ( A . 1 4)  into 
Eqs.  (A. 9)  and  (A. 10)  and  comparing  Eqs.  (A. 5)  and  (A. 8)  with  Eqs.  (A.1) 
and  ( A.4 ) ,  respectively,  we  obtain,  for  below  the  cloud, 


F+(z)  -  Fb(z)  ,  (A. 15) 

F+(z)  -  [rC  F+ ( zfc)  +  tG  F+(zt)  -  (l-eC)  otJ]  [l  -  ef(zb-z,Tb)]  +  F*(z), 

(A. 16) 


and  for  above  the  cloud, 


F+(z)  =  [rC  F+(zt)  ♦  tC  F+(zb) 
F+(z)  =  F*(z)  . 


(l-EC)  oTb]  [1  -  ef(z-zt,Tt)]  +  Fb(z)  , 

(A. 17) 
(A. 18) 


The  above  modification  replaces  complicated  and  redundant  computations  of 
the  transmitted  and  reflected  fluxes  with  one-time  calculations  of  input 
fluxes  at  the  cloud  boundary,  and  thus  reduces  the  computation  time  by 
one  half.  The  resulting  radiative  heating  rate  and  flux  profiles  from 
this  modified  scheme  differ  from  those  in  the  original  scheme  by  a 


negligible  amount. 


Appendix  B 


DOCUMENTATION  OF  RADIATION  PARAMETERIZATION 
PROGRAMS 

B. 1  General  Description  of  the  Program 

The  basic  principles  for  calculating  both  IR  and  solar  radiation 
fluxes  and  heating  rates  have  been  described  in  Sections  2  and  3.  Here, 
we  document  detailed  computer  routines  for  IR  and  solar  parameterization 
programs  for  use  In  general  circulation  models. 

The  entire  radiation  program  is  composed  of  three  subroutines 
denoted  as  PRERAD,  IRRAD,  and  IRCLR.  The  subroutine  PRERAD  serves  as  a 
driver  program  that  connects  the  main  part  of  radiation  calculations  with 
other  GCM  programs.  It  Interpolates  the  model-generated  temperature  and 
humidity  profiles  from  the  main  part  of  the  GCM  on  layers  designed  for 
radiation  calculations .  The  cloud  position,  along  with  cloud  cover  and 
thickness,  are  determined  and  used  as  input  to  the  subroutine  IRRAD, 
which  is  subsequently  called.  The  computed  radiative  heating  rates  for 
radiation  layers  from  the  subroutine  IRRAD  are  then  weighted  by  density 
to  get  the  mean  radiative  heating  rate  for  each  GCM  layer. 

The  subroutine  IRRAD  is  the  main  program  for  radiation  calculations. 
It  is  composed  of  IR  and  solar  parts.  In  the  IR  part,  the  path  lengths 
in  pressure  coordinates  are  first  computed  and  then  adjusted  in  accor¬ 
dance  with  a  number  of  temperature  and  pressure  correction  schemes 
denoted  in  Section  2.  Next,  the  coefficients  of  the  third-order  poly¬ 
nomial  for  IR  emissivities  are  obtained  by  interpolation  between  the 


prescribed  reference  temperatures  listed  in  Table  1 .  The  IR  program  is 
then  connected  to  a  large  computational  loop  in  which  radiative  fluxes 
for  each  cloud  configuration  are  calculated.  The  strategy  for  computing 
fluxes  has  been  illustrated  in  Fig.  3.  For  computing  fluxes  in  a  clear 
atmosphere,  the  subroutine  IRCLR  is  called  with  the  lower  and  upper 
boundaries  being  the  surface  and  top  of  the  atmosphere,  respectively. 

For  computing  fluxes  above  and  below  a  high  cloud  deck,  we  first 
obtain  the  flux  profiles  for  both  above  and  below  the  cloud  by  calling 
the  subroutine  IRCLR,  assuming  that  the  cloud  is  a  blackbody.  Then  the 
fluxes  are  adjusted  to  account  for  the  effect  of  transmission  through  and 
reflection  by  the  high  cloud,  according  to  Eqs.  (A. 1 0)- (A. 1 3) .  On  the 
other  hand,  fluxes  above  a  black  cloud  without  an  overhead  high  cloud,  or 
between  two  black  clouds,  or  below  a  black  cloud  are  obtained  by  calling 
IRCLR.  Intracloud  flux  profiles  are  obtained  either  by  interpolation 
(high  cloud)  or  by  the  blackbody  assumption  (middle  and  low  clouds). 

The  solar  program  can  be  divided  into  three  parts.  First,  the  path 
lengths  of  each  model  layer  for  f^O,  CC^,  and  0^  are  computed.  These 
layered  path  lengths  are  then  substituted  into  equations  for  the  computa¬ 
tion  of  gaseous  absorpti vlties .  Second,  fluxes  for  the  clear  sky  are 
obtained.  Third,  the  computation  of  cloudy  sky  fluxes  is  further  divided 
into  three  parts,  according  to  the  number  of  cloud  decks  Involved.  The 
first  part  consists  of  computing  the  cloud  optical  properties  based  on 
the  input  cloud  thickness  and  prescribed  cloud  water  content.  In  the 
event  that  there  is  only  one  cloud  deck,  the  flux  computation  is  carried 
out  based  on  the  ray  tracing  principle.  The  other  two  parts  consist  of 
computing  two-  and  three-deck  fluxes,  according  to  the  methods  described 


In  Section  3. 


Finally,  the  subroutine  IRCLR  computes  IR  fluxes  In  clear  regions. 
It  can  be  used  to  compute  the  entire  clear  column  between  the  model  top 
and  surface,  the  clear  region  between  the  low  and  middle  clouds,  as  well 
as  between  these  clouds  and  the  surface.  To  carry  out  the  flux  integra¬ 
tion  for  clear  columns  the  flux  equation  is  generalized  in  the  following 
f lnlte-dlfference  form 


elm,J  +  ^  ^J+1/2  ^elm,J 

J"ns 


-wM  •  <B-’> 


B.2  Interfacing  with  the  AFGL  GCM 

A  schematic  flow  diagram  illustrating  the  connection  between  the 
AFGL  GCM  and  the  radiation  program  is  depicted  in  Fig.  B.l.  In  the  AFGL 
GCM  main  program  the  subroutine  LALOOP  is  called  to  compute  large-soale 
parameters  for  a  specific  time  step.  The  subroutine  NLPROD  is  called  in 


Table  B.1  Assignments  of  parameters  involved  in  Eq.  (B.1). 


1  mb  +  1  m  mb  -1  oTb 

2  m  +  1  mb  mt  1  oTb 

mb:  level  index  for  the  base  of  the  integration  domain 

mb:  level  index  for  the  top  of  the  integration  domain 

Tb:  temperature  at  the  base  of  the  integration  domain 

Tfc:  temperature  at  the  top  of  the  integration  domain 


LALOOP  to  compute  the  nonlinear  product  In  LALOOP.  In  NLPROD,  the  sub¬ 
routine  SUBRAD  is  called  to  compute  radiative  transfer  parameters.  In 
the  subroutine  SUBRAD,  the  solar  zenith  angles  are  computed  through  the 
subroutine  ZENITH,  the  cloud  cover  field  is  computed  through  the  sub¬ 
routine  CLOUD,  and  then  the  subroutine  PRERAD  is  called  to  compute  the 
radiative  heating  rate,  which  is  subsequently  added  to  the  large-scale 
heating  rate  previously  computed  in  NLPROD.  The  vertical  grid  structure 
used  in  radiation  routines  is  consistent  with  the  AFGL  GCM  grid  struc¬ 
ture,  so  that  there  is  no  vertical  interpolation  of  the  large-scale 
parameters  required  In  the  radiation  routine  and  computation  time  is 
saved.  The  effects  of  the  vertical  resolution  have  been  shown  in  Section 
^  through  the  Intercomparison  of  clear  and  cloudy  heating  rate  profiles. 

Another  link  between  the  AFGL  GCM  and  radiative  routines  Is  the 
computation  of  the  surface  radiative  flux  balance.  The  downward  IR  flux 
and  downward  net  solar  flux  at  the  surface  for  each  grid  are  generated  in 
radiative  routines  and  transferred  to  the  AFGL  GCM  through  the  common 


block  VARIA ,  for  computing  the  surface  energy  balance  In  the  routines  for 


LALOOP 


i 


I 


g.  B.1  A  schematic  flow  diagram  Illustrating  the  Interface  of  the 
AFGL  GCM  with  the  current  radiation  programs.  The  routines 
contained  in  the  dashed  lines  are  radiation  routines. 


boundary  layer  fluxes. 


B.3  Changes  Required  in  the  Radiative  Routines  for  Converting  Variable 


Resolution  to  Fixed  Resolution 


The  radiative  routines  have  been  modified  to  accomodate  computations 


using  variable  vertical  resolutions,  since  a  variable-resolution  version 


of  the  AFGL  GCM  is  already  available.  Every  time  the  variable-resolution 


radiative  routines  are  converted  to  the  fixed-resolution  version,  certain 


changes  must  be  made.  These  changes  can  be  divided  into  three  groups  as 


follows: 


1)  The  common  blocks  transfer  parameters  between  the  AFGL  GCM  and 


radiative  routines:  these  common  blocks  include  CCONST,  GRDPT,  ALEG, 


SURFAC ,  VRTST,  BL,  CZALB,  STORE  and  BLFXT.  The  specif ications  of  array 


dimensions  need  to  be  changed  so  as  to  be  consistent  with  the  AFGL  GCM. 


2)  The  DATA  statements  assign  values  to  JUDGE  arrays:  JUDGE  (J,K) 


is  a  2-dimenslonal  array.  It  is  introduced  to  the  program  for  prescrib¬ 


ing  cloud  top  and  base  level  indices  and  cloud  types.  The  level  Indices 


are  usej  in  radiative  routines  to  specify  clear  regions  separated  by 


:louds.  The  cloud  types  are  required  information  for  computing  the  solar 


radiative  properties  of  clouds.  The  first  argument  J  indicates  the  cloud 


configuration.  The  high,  middle,  and  low  cloud  decks  are  given  indices 


of  1,  2,  and  3,  respectively,  whenever  that  cloud  deck  is  present  in  the 


grid.  If  no  cloud  deck  is  present,  the  index  is  zero.  Then  the  value  of 


J  is  obtained  by 


l  [I(j)]‘ 
J-1 


(B.3) 


I 


where  I(j)  is  the  index  for  the  jth  cloud  deck.  The  second  argument  K 
denotes  the  cloud  top  (K  =  1,4,7)  and  base  (K  -  2,5,8)  level  indices,  and 
the  cloud  type  (K«  3,6,9).  For  a  two-deck  (high  and  low  clouds)  case, 

J  =  1  +  9  =  10 

JUDGE  (10,1)  is  the  high  cloud  top  level  index,  JUDGE  (10,2)  the  high 
cloud  base  level  index,  JUDGE  (10,3)  the  high  cloud  type  index,  JUDGE 
(10,1))  the  low  cloud  top  level  index,  JUDGE  (10,5)  the  low  cloud  base 
level  index,  JUDGE  (10,6)  the  low  cloud  type  index,  and  JUDGE  (10,7)  - 
JUDGE  (10,8)  =  JUDGE  (10,9)  -  0.  These  quantities  of  JUDGE  need  to  be 
respecifled  every  time  the  vertical  resolutions  are  changed. 

3)  The  DATA  statements  assign  humidity  and  ozone  mixing  ratios: 
the  variable  El  is  the  humidity  mixing  ratio.  The  AFGL  GCM  only  provides 
humidity  values  up  to  a  certain  sigma  value  (usually  corresponding  to  the 
tropopause).  Thus  the  humidity  mixing  ratio  in  the  stratosphere  needs  to 
be  respecified,  according  to  climatology,  when  the  vertical  resolution  is 
changed.  The  variable  Q1  is  the  ozone  mixing  ratio.  Since  the  AFGL  GCM 
does  not  provide  the  ozone  mixing  ratio,  this  quantity  also  needs  to  be 
prescribed  over  the  whole  vertical  column  according  to  the  climatological 
mean.  Finally,  the  quantity  LT0P,  which  denotes  the  top  level  index  for 
the  Rayleigh  reflecting  layer,  needs  to  be  changed. 


Appendix  C 


DIRECTORY  OF  ARRAYS  AND  VARIABLES  IN 
THE  RADIATIVE  ROUTINES 


Subroutine  PRERAD 

Atmospheric  properties 

Array 

CL9  (8) 

El  (HO) 

Q1  (HO) 

T1  (HO) 

PTX  (6) 

W  (96,2,12) 

T  (96,2,12) 

SL  (12) 

TO  (12) 

Variable 

HBAR 

HQo 

Storage  arrays 
Array 

FNIR  (96,76) 

FNSOL  (96,76) 

KFLUX  (96,76,9) 

RIIR  (96,2,13) 

RISOL  (96,2,13) 

IFLUX  (76,9) 


Meaning 

sigma  structure  array 
layer  humidity  mixing  ratio 
layer  ozone  mixing  ratio 
layer  temperature 
surface  pressure  and  temperature 
humidity  field  (longitude  •  sigma)  for 
two  equator-symmetric  latitude  circles 
from  the  GCM 

temperature  field  (longitude  •  sigma)  for 
two  equator-symmetric  latitude  circles 
from  the  GCM 

layer  sigma  values  from  the  GCM 
reference  temperature 

Meaning 

prescribed  relative  humidity  in  the 
stratosphere 

saturated  humidity  mixing  ratio  in  the 
stratosphere 


Meaning 

net  upward  IR  flux  field  at  the  surface 
net  downward  solar  flux  field  at  the 
surface 

storage  array  for  IR,  solar,  and  net 
fluxes  at  the  top  of  the  atmosphere, 
within  the  atmosphere  and  the  surface 
IR  radiative  heating  rate  for  two 
equator-symmetric  latitude  circles 
solar  radiative  heating  rate  for  two 
equator-symmetric  latitude  circles 
storage  array  for  zonally  averaged  IR, 
solar,  and  net  flux  fields 


3.  Radiative  parameters 


Array 

H9  (4) 

H  (40) 

RIRAD  (96,2,13) 


SZ  (96,2) 

SALB  (96,2) 

Variable 

UIRF 

DIRF 

CZ 

SA 

DSOLF 

4.  Cloud  parameters 

Array 

CB  (7,3) 

CT  (7,3) 

FRACT  (7) 

ICL  (96,76,4) 

L23  (6.7) 

JUDGE  (14,9) 


LTYPE  (7,3) 

Variable 

LCBOT 

LCTOP 

IFRACT 

ISUM 

5.  Constants 

Array 

PHI  (76) 
LNGTD  (96) 
GLAT  (38) 


Meaning 

solar  heating  rate 
IR  heating  rate 

net  radiative  heating  (longitude  •  sigma) 
for  two  equator-symmetric  latitude 
circles 

cosine  of  the  solar  zenith  angle  field 
surface  albedo  field 

Meaning 

net  upward  IR  flux  at  the  surface 
downward  IR  flux  at  the  surface 
cosine  of  the  solar  zenith  angle  for  one 
grid  point 

surface  albedo  for  one  grid  point 
downward  solar  flux  at  the  surface 


Meaning 

cloud  deck  base  level  sigma 
cloud  deck  top  level  sigma 
cloud  cover 

cloud  cover  field  In  decades  of 
percentage 

cloud  deck  top  and  base  level  Indices 
reference  array  storing  cloud  deck  top 
and  base  level  indices  and  cloud  type 
index  for  all  configurations 
cloud  type  index 

Meaning 

high  cloud  base  level  Index 

high  cloud  top  level  index 

cloud  cover  for  one  grid  point  in  decades 

of  percentage 

sum  of  square  of  cloud  deck  indices 


Meaning 

grid  point  latitude  value 
grid  point  longitude  value 
Gaussian  latitude  values 
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Constants  (cont.) 


Variable 

SECPDY 

SIG 

PII 

NCLD 

NCLDM 

NGCM 

MP 

NP 

NP2 

NP2P 

KP 

KPHL 

KPM1 

KPP1 

LTOP 


Meaning 

seconds  per  day 

Stef an-Boltzman  constant 

TT 

total  number  of  cloud  configurations 
(Including  clear  case)  In  a  grid 
NCLD  -  1 

total  number  of  GCM  model  levels 
total  number  of  GCM  longitudinal  grids 
total  number  of  GCM  latitudinal  grids 
NP/2 
NP2  +  1 

total  number  of  GCM  vertical  grids 
total  number  of  humidity  data  void  layers 
+  1 

KP  -  1 
KP  +  1 

level  Index  for  the  top  of  the  Rayleigh 
scattering  layer 


Subroutine  IRRAD 

Atmospheric  properties 

Arrays 

PRM  (40  ) 

TEM  (40) 

AVGQ  (40) 

F  v 5 ,  40) 

URF  (5,  40) 

UWCOM  (6,  40) 

PTX  (6) 

T1  (40) 

El  (40) 

Q1  (40) 

CL9  (80) 

PBL  (40) 

PBY  (40) 

PBD  (40) 

PMD  (40) 

PMU  (40) 

VPY  (40) 

0PY  (40) 

CPY  (40) 

VDS  (40) 

ODS  (40) 


Meaning 

model  level  pressure 
model  level  temperature 
model  level  humidity 

pressure  and  temperature  scaling  factor 
pressure  and  temperature  scaled  path 
unsealed  path 

surface  pressure  and  temperature 
layer  temperature  from  MAIN 
layer  humidity  mixirg  ratio  from  MAIN 
layer  ozone  mixing  ratio  from  MAIN 
sigma  structure  array 
model  level  pressure 
model  layer  pressure 
model  layer  pressure  thickness 
mean  scaled  pressure  (counted  +) 
mean  scaled  pressure  (counted  +) 
model  layer  H20  path  thickness 
model  layer  0,  path  thickness 
model  layer  C02  path  thickness 
Integrated  path  for  H20  (downward) 
integrated  path  for  03  (downward) 


1.  Atmospheric  properties  (cont.) 


Arrays 


Meaning 


CDS  (40) 
VUS  (40) 
OUS  (40) 
CUS  (40) 


Integrated  path  for  C02  (downward) 

integrated  path  for  H20  (upward) 

integrated  path  for  03  (upward) 

integrated  path  for  C02  (upward) 


Variables 


Meaning 

/  uw  dpw,  pw:  water  path 
/  dpw 


2.  IR  gaseous  radiative  properties 


Arrays 

ABSB  (40,2,6) 

REFT  (4) 

RK  (5,40) 

RJ  (5,40,4) 

EE  (5,4,4) 

INDEXT  (40) 

DT  (40) 

CDT  (40) 


Meaning 

emissivity  between  the  surface  and  level 
of  interest 

reference  temperature  for  determining 
emissivity 

reference  emissivity  for  path  lengths 
less  than  10-7  g  cm"* 
temperature- Interpolated  emissivity 
coefficients 

broadband  emissivity  coefficients  for 
five  bands  and  four  temperatures 
reference  temperature  index  for  each 
model  level 

interpolation  factor  for  temperature 
1  -  DT 


Variable 


Meaning 
INDEX (M )  +  1 


3.  Solar  gaseous  radiative  properties 


Arrays 

Meaning 

WVA 

(40) 

/  Avdv  for  downward  direct  solar  radia¬ 
tion  in  2.7  urn  band 

WVB 

(40) 

/  Avdv  for  downward  diffuse  sola’"  radia¬ 
tion  in  2.7  um  band 

WVC 

(40) 

/  A  dv  for  upward  diffuse  solar  radiation 
in  2.7  ym  band 

WEI 

(40) 

weighting  factor  e  for  downward  direct 
solar  radiation  in  2.7  ym  bend 

WE  J 

(40) 

e  for  downward  diffuse  solar  radiation  in 
2.7  ym  band 
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Solar  gaseous  radiative  properties  (cont.) 


Arrays 

WEK  (40) 

BLS  (7) 

CLS  (7) 

DLS  (7) 

ELS  (7) 

FLS  (7) 

FOZ  (18) 
AOZ  (18) 
ACD  (40) 
ACY  (40) 
ACU  (40) 
AVD  (6,40) 
AVY  (6,40) 
AVU  (6,40) 
ABD  (40) 

ABY  (40) 

ABU  (40) 

AAD  (40) 

AAY  (40) 

AAU  (40) 

AOU  (40) 

AOx  (40) 

AOD  (40) 

AUX  (6) 

AUX  (6) 

Variables 


Meaning 

e  for  upward  diffuse  solar  radiation  in 
2.7  pm  band 
Av 

C-factor 

D-factor 

K/D-factor 

fraction  of  H20  bands 
fraction  of  03  bands 
03  absorption  coefficients 
C02  absorption  (downward  direct) 

C02  absorption  (downward  diffuse) 

C02  absorption  (upward  diffuse) 

H20  absorption  (downward  direct) 

H20  absorption  (downward  diffuse) 

H20  absorption  (upward  diffuse) 
broadband  absorption  for  H20  +  03  (down¬ 
ward  direct) 

broadband  absorption  for  H20  +  03  (down¬ 
ward  diffuse) 

broadband  absorption  for  H20  +  03  (upward 
diffuse) 

broadband  absorption  for  H20  (downward 
direct) 

broadband  absorption  for  H20  (downward 
diffuse) 

broadband  absorption  for  H20  (upward 
diffuse) 

broadband  absorption  for  03  (upward 
diffuse) 

broadband  absorption  for  03  (downward 
diffuse ) 

broadband  absorption  for  03  (downward 
direct ) 

constants  for  converting  scaled  path  to 
logarithmic  path 

Meaning 

10C/D 

1-A+(zti ) 

1-A+ (zt2) 

1-A+(zto) 

1-A*(zti ) 

1  — A ♦ ( ^ 

1  -  A+  ( z^  o ) 

1-A+(zt1 ) 

1-A+(zt?) 

1 -A+ ( zt3 ) 


(upward 

(downward 


(downward 


IR  cloud  radiative  properties 


Arrays 


Meaning 


YRE  (6) 
YTR  (6) 


YEM  (6) 


Variables 


coefficients  for  high  cloud  IR  reflec¬ 
tivity 

coefficients  for  high  cloud  IR  transmis¬ 
sivity 

coefficients  for  high  cloud  IR  emisslvlty 
Meaning 


high  cloud  IR  emlssivity 
high  cloud  IR  reflectivity 
high  cloud  IR  transmissivity 


Solar  cloud  radiative  properties 


Arrays 
Y  (8,4*1) 

YT  (8,3,3) 
YA  (4,4,4) 
TREF  (4) 


RC  (2) 
TC  (2) 


RCBAR  (2) 


TCBAR  (2) 


ACBAR  (2) 
AC  (2) 


WREF  (4) 


Meaning 

coefficients  for  solar  cloud  radiative 
properties 

coefficients  for  solar  cloud  radiative 
property  with  temperature  adjustments 
coefficients  for  solar  cloud  absorption 
properties 

reference  temperature  for  computing  tem¬ 
perature  adjustment  for  each  type  of 
cloud 

direct  solar  cloud  reflectivity  for  1-  or 
2-layer  clouds 

direct  solar  cloud  transmissivity  for  1- 
or  2-layer  clouds 

diffuse  solar  cloud  reflectivity  for  1- 
or  2-layer  clouds 

diffuse  solar  cloud  transmlssl vlty  for  1- 
or  2-layer  clouds 

diffuse  solar  cloud  absorptivity  for  1- 
or  2-layer  clouds 

direct  solar  cloud  absorptivity  for  1-  or 
2-layer  clouds 

maximum  cloud  liquid  water  content 


Variables 


Meaning 


TTTT 

RRRR 

AAAA 

AAA 

RRR 

TTi 


cloud  transmittance 
cloud  reflectance 
cloud  absorption 
absorptivity  a 
reflectivity  r 
transmissivity  t 


5.  Solar  cloud  radiative  properties  (cont.) 


Variables 


Meaning 


RCTC 

RCL 

TCL 

RUO 

TUO 

RCBR1 

RCBR2 

RCL1 

RCL2 

RCTC1 

RCTC2 

RTBR 

TCBR1 

TCBR2 

TCL1 

TCL2 


r  +  t  +  a 
normalized  r 
normalized  t 
normalized  r 
normalized  t 
RCBARO  )/RCTC1 
RCBAR (2)/RCTC2 
RC(1 )/RCTC1 
RC(2)/RCTC2 
r?  +  t?  +  a? 
r?  +  tS  +  af 
?k  +  tg  ♦  a£ 
TCBARO  )/RCTC1 
TCBAR(2)/RCTC2 
TC(1  )  /RCTC1 
TC(2 )/RCTC1 


6.  IR  radiative  transfer  quantities 


Arrays 


Meaning 


ST  (8,40) 

NX  (4) 

NCHAN  (8,3) 

AT  (40) 

SS  (6,40) 

TQD  (8) 

SGN  (8) 

AUX  (6) 

AUY  (6) 

H  (80,40) 


components  of  radiative  flux  for  various 
cloud  types  plus  weighted  average 
index  for  components  above  high  cloud 
flag  matrix  for  variable  INDD  elements 
(1  =  variable) 

T1  **  (layer  temperature) 

Components  of  radiative  fluxes  for  one 
band 

Tg  (Tg :  surface  temperature) 
prefix  matrix 

constants  for  converting  scaled  path  to 
logarithmic  path 
IR  cooling  rate  vector 


Variables 


Meaning 


W1 

W2 

DFLUX 


linear  interpolation  factor  for  in-cloud 

radiative  fluxes 

(1-cc)  [F+(zt)-oT&)  +  rc  F+(zb) 

[Eqs .  (A. 16),  (A. 17)] 

downward  IR  flux  at  the  surface 

net  upward  IR  flux  at  the  surface 


7.  Solar  radiative  transfer  quantities 
Arrays  Meaning 


H9  (40) 
FDS  (40) 
FUS  (40) 

Variables 

CZ 

SA 

ANGLE 

SOLAR 

ALB 

EINSO 

RAS 

RAL 

REF 

FZT 

FZBU 

T 

TD 

TU 

SI  2 

S23 

S34 

AXX 

AYY 

DSOLF 

G1D 

G2D 

G3D 

G4D 

G1U 

G2U 

G3U 

G4U 

G1Y 

G2Y 

G3Y 

G4Y 

XG1 

XG2 

YG1 

YG2 

YG3 


solar  heating  rate 
downward  solar  flux 
upward  solar  flux 

Meaning 

solar  zenith  angle  y0 
solar  surface  albedo  rg 
solar  zenith  angle  p0 
solar  constant  S0 
solar  surface  albedo  r„ 

S0Mo 

Rayleigh  layer  reflectance  (diffuse) 
Rayleigh  layer  reflectance  (direct) 
total  atmospheric  Rayleigh  layer  reflec¬ 
tance 

downward  flux  at  the  top  of  Rayleigh 
layer 

upward  flux  at  the  cloud  base 
1  -  [ABD  (layer)  -  ABD  (LBUP)]/[1  -  ABD 
(LBUP) ] 

1  -  [ABY  (layer) 

( LBUP ) ] 

1  -  ABU  (LBCL) 

Il2 

|23 
~34 

1  - 


-  ABY  (LBUP ) ]/[ 1  -  ABD 


ABU 

ABD 


( LTCL) 
(LBUP) 


downward  net  solar  flux  at  the  surface 

G1  + 

G2* 

Go* 

oj* 

2r 

G2t 

h* 

G^t 

Sl  + 

Go  + 

xi 


yi 

y  2 

y3 


s 

V 

a 


* 

ft 

M 

$ 


5 

s 

X 


V  j 

vl 

* 

*  r 

vl 

i 

» 

5  * 
rJ 

s 

.V 


I 


a 

s' 

y 


81 


ppnonCTTT 

s 


> 

j 


7.  Solar  radiative  transfer  quantities  (cont.) 
Variables  Meaning 


i 


YGM 

yn 

TOD 

V 

T1D 

V 

V 

T2D 

T3D 

h* 

T1U 

h* 

T2U 

I2* 

T3U 

l3  + 

T1Y 

Il  + 

T2Y 

l2  + 

T3Y 

Cloud  properties 

V 

j 

; 

i 


3 


i 


Arrays 


LC  (6) 

NLC  (6) 

L23  (6,7) 
FRACT  (7) 

THK 

WTCNT  (6,7,3) 
LTYPE  (7,3) 
Variables 


ISUM 

LCTOP 

LCBOT 

WZ 

PART 

LTCL 

LBCL 

DZ 

IC 

ACT 

ICX 

ICY 

LTUP 

LBDN 

LTDN 

LBUP 

LTUP1 


Meaning 

Cloud  top  and  cloud  base  level  indices 
(counted  top  ♦) 

cloud  top  and  cloud  base  level  indices 
(counted  bottom  +) 

cloud  top  and  base  level  index  (top  *) 
cloud  fraction 
cloud  thickness 

liquid  water  content  for  different  types 

of  cloud3 

cloud  type  index 

Meaning 


summation  of  square  of  cloud  deck  index, 

first  argument  of  JUDGE 

cirrus  cloud-top  index  (+) 

cirrus  cloud-base  index  (+) 

total  liquid  water  path 

cloud  fraction 

cloud  top  index 

cloud  base  index 

cloud  thickness 

cloud  type  index 

average  temperature  in  cloud 

(cloud  type  index)  *  2-1 

(cloud  type  index)  *  2 

LTCL  -  1 

LBCL  +  1 

LTCL  +  1 

LBCL  -  1 

LC(1 )  -  1 

82 


l 

1 

a 

x 


a 


^1 

w 

n 

l 


% 


\ 


A 


A 


73 

s 

f 

ft 
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8.  Cloud  properties  (cont.) 


Variables 

Meaning 

LTUP2 

LC( 3)  -  1 

LTUP3 

LC(5)  "  1 

LBDN1 

LC(2)  +  1 

LBDN2 

LC(4)  +  1 

LBDN3 

LC ( 6 )  +  1 

LTCL1 

LC(1  ) 

LTCL2 

LC  ( 3 ) 

LTCL3 

LC(5 ) 

LBCL1 

LC(2) 

LBCL2 

LC(4) 

LBCL3 

LC(6) 

LBUP1 

LC(2)  -  1 

LBUP2 

LC(4)  -  1 

LBUP3 

LC(6)  -  1 

LTDN1 

LC ( 1 )  +  1 

LTDN2 

LC  ( 3 )  +  1 

LTDN3 

LC(5)  +  1 

NMM 

NLC(1 )  -  1 

NPP 

NLC ( 2 )  +  1 

Constants 

Variables 

Meaning 

LEVEL 

total  number  of  model  levels 

NLYR 

total  number  of  model  layers 

NCLD 

total  number  of  cloud  types 

LTOP 

index  for  the  top  level  of  Rayleigh's 
reflection  layer 

SIG 

Stefan-Boltzmann  constant 

CPR 

constant  for  computing  reduced  C02  path 

NLVL 

total  number  of  model  levels 

NTOT 

total  number  of  non-overlap  bands 

NTOTAL 

total  number  of  bands 

LSRFC 

surface  or  cloud  top  level  index  below 
the  high  cloud 

NCLDP 

total  cloud  type  number  +  1 

LAYER 

total  number  of  model  layers 

LAYERM 

LAYER  -  1 

CP 

specific  heat  at  constant  pressure  for 
air 

IFLG 

flag  for  printing  radiative  transfer 
quantities 

KKE 

number  of  cloud  decks 

C.  Subroutine  IRCLR 


1 .  Atmospheric  variables 


Arrays 


UPATH  (6,40) 


URF  (5,40) 

AT  (40) 

UWCOM  (6,40) 
TEM  (40) 


2.  Radiative  Variables 


Arrays 


ABSB  (40,2,6) 


EMMI  (6,40) 
FLUXI  (6,6,40) 
SS  (6,40) 

DRK  (4) 

DEMIS  (6) 

RK  (5,40) 


RJ  (5.40,4) 


INDD  (8,4) 
TQD  (8) 

NCHAN  (8,  3) 


SCM  (8) 
AUX  (6) 
AUY  (6) 


3.  Constants 


Variables 


SIG 

CPR 

NLVL 

NTOT 

NT0TAL 

LEVB 

LEVT 

NNS 

NNE 

NN3 

NN4 

NNSP 


Meaning 


path  length  between  emission  and  contri¬ 
buting  levels 
path  length 

T1 "  (layer  temperature) 
unadjusted  path  length 
model  level  temperature 


Meaning 


emlsstvlty  between  surface  and  level  of 
Interest 

emlsstvlty  corresponding  to  UPATH 

components  of  fluxes 

components  of  fluxes 

a,  u1"1  (emlsstvlty  term)  r 

1  -  Ef (u,Ts)  or  e”(u+Au,T)  -  el(u,T) 

reference  emlsslvlty  for  path  lengths 

less  than  10‘7  g  cnr2 

temperature-interpolated  emlsslvlty  coef¬ 
ficients 

prefix  and  lndlce  matrix 

Tg  (Tg:  surface  temperature) 

flag  matrix  for  variable  INDD  elements 

1  -  variable) 

prefix  matrix 

constants  for  converting  scaled  path  to 
logarithmic  path 


Meaning 


o:  Stefan-Boltzmann  constants 
constants  for  computing  C02  path 
total  number  of  model  levels 
total  number  of  non-overlap  bands 
total  number  of  bands 
level  Index  of  clear  region  base 
level  Index  of  clear  region  top 
emission  level  index 
contributing  level  index 
some  reference  index 
boundary  index 
NNS  +  1 


Si 


8 


I 


kS 


» 

ft 


i 


% 

w 
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